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TECHNICAL NOTE 3096 

THEORETICAL PREDICTION OF PRESSURE DISTRIBUTIONS ON 
NONLIFTING AIRFOILS AT HIGH SUBSONIC SPEEDS 
By John R. Spreiter and Alherta Alksne 


SUMMARY 


Theoretical pressure distributions on nonlifting circular-arc air- 
foils in two-dimensional flows with high subsonic free -stream velocity 
are found by determining approximate solutions , through an iteration 
process, of an integral equation for transonic flow proposed by 
Oswatitsch. The integral equation stems directly from the small- 
disturbance theory for transonic flow. This method of analysis pos- 
sesses the advantage of remaining in the physical, rather than the 
holograph, variables and can be applied to airfoils having curved sur- 
faces. After discussion of the derivation of the integral equation and 
qualitative aspects of the solution, results of calculations carried 
out for circular-arc airfoils in flows with free -stream Mach numbers up 
to unity are described. These results indicate most of the principal 
phenomena observed in experimental studies. At subcritical Mach numbers, 
the pressure distribution is symmetrical about the midchord position 
and the drag is zero. The magnitude of the pressure coefficient is 
found to increase more rapidly with increasing Mach number than the 
Prandtl-Glauert rule would indicate.. When the critical Mach number is 
exceeded, compression shocks occur, the fore-and-aft symmetry of the 
pressure distribution is lost, and the airfoil experiences a drag force. 
As the Mach number is increased further, the shock wave becomes of 
greater intensity and moves rearward along the chord, thereby producing 
a rapid increase in the magnitude of the pressure drag coefficient. At 
Mach numbers close to unity, the variation of the pressure, local Mach 
number, and drag conforms, within the limitations of transonic small 
perturbation theory, to the known trends associated with the Mach number 
freeze. Some comparisons with experimental results are also included. 

The solutions are obtained using an iteration process which differs 
from the classical methods in that the quadratic nature of the integral 
equation is recognized. If the iteration calculations are started using 
the linear-theory solution, it is shown that the retention of the quad- 
ratic feature has the interesting effect of forbidding shock-free super- 
critical second-order solutions. In order to obtain solutions for 
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supercritical Mach numbers, it is necessary to start the iteration cal- 
culations with a velocity or pressure distribution which contains a 
compression shock. When this is done, it is found that the iteration 
procedure converges to a definite result. 


INTRODUCTION 


The theoretical problem of transonic flow about thin wings has 
been discussed by numerous authors in recent years. Since the basic 
equations are nonlinear and of mixed type, the difficulties are great, 
and progress has been made only through expenditure of considerable 
effort. At first, only the basic equations and the similarity rules 
were established. (See refs. 1 through 12.) More recently, a small 
number of actual solutions have been determined. At the present time, 
the most complete theoretical results are those of Guderley and Yoshihara, 
Vincenti and Wagoner, Cole, and Trilling (refs. 13 through l8) for the 
flow about wedge airfoils at both subsonic and supersonic speeds. These 
were all obtained by transforming the equations to hodograph variables 
whereby the differential equation becomes linear although still of mixed 
type. Superposition of solutions is then possible, but the boundary 
conditions generally become very complicated. It is because of the 
latter difficulty that all the solutions mentioned above are for wedge 
sections. A further disadvantage of the hodograph method is that it is 
definitely restricted to two-dimensional flows, there being no known 
transformation which linearizes the equation for three-dimensional com- 
pressible flows. 

If the hodograph transformation is not introduced, there are avail- 
able no direct methods of solution. However, various iteration methods 
have been used to study flows with high subsonic free-stream velocities. 
(See ref. 19 for a resume.) Almost all these account for the compressi- 
bility effects by source distributions throughout the flow field and 
start with either the solution for incompressible flow or for linearized 
compressible flow as the first approximation. A second approximation 
is calculated from the first and so on. It was not found possible, how- 
ever, to iterate starting with a typical shock-free subsonic flow solu- 
tion and obtain a typical transonic flow field in which the supersonic 
region ends with a shock. 

Oswatitsch has presented another method in references 20 and 21 for 
determining the transonic pressure distribution on thin airfoils in flows 
with subsonic free-stream velocity. The analysis is carried out in the 
physical rather than the hodograph variables, and leads to a nonlinear 
integral equation in which the unknown velocity appears outside as well 
as inside the integral. Oswatitsch finds approximate solutions not by 
iteration, but by introducing various functions containing undetermined 
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parameters into the integral equation and by determining the parameters 
so that the integral equation is satisfied at a small number of points 
on the airfoil. The method is applied to determine the pressure distri- 
bution on circular -arc and NACA symmetrical four-digit airfoils. The 
results show certain characteristics of transonic flow such as the 
appearance of shock waves and their rearward movement across the chord 
with increasing Mach number. It is disconcerting, however, that the 
method fails to give proper results at high subsonic Mach numbers 
(greater than about 0.88 for 6-percent-thick circular-arc sections), 
provides a multiplicity of solutions at supercritical Mach numbers, and 
permits the integral equation to be satisfied at only a very limited 
number of points. 

The present work is based on the integral equation of Oswatitsch 
but an iteration process is used to obtain approximate solutions. This 
procedure permits the integral equation to be satisfied at a much larger 
number of points than in the original method of Oswatitsch, gives approxi 
mate solutions at all Mach numbers up to unity, and appears to avoid any 
multiplicity of solutions. The method is applied to determine the theo- 
retical pressure distribution on symmetrical circular-arc airfoils at 
zero angle of attack. Except for phenomena that are primarily of viscous 
origin, such as boundary-layer separations, etc., these results exhibit 
most of the experimentally observed features of transonic flows. 

Attention is also called to reference 22 by Gullstrand, an associate 
of Oswatitsch, in which transonic flows about thin airfoils are investi- 
gated by still another extension of Oswatitsch 1 s integral-equation theory 
Gullstrand sought to determine approximate solutions by iteration, 
although his procedure differs considerably from that described herein. 
His method, however, succeeded in determining solutions only when the 
Mach number was less than about 0.90 for 6 -percent -thick sections. 

Results were given for the pressure distribution at sonic speed in a 
second paper by Gullstrand (ref. 23 ), but they were obtained by introduc- 
ing a new and more complicated integral equation than that of Oswatitsch 
used herein. In contrast to the present analysis in which the entire 
solution is obtained from the integral equation, Gullstrand uses the 
integral equation to determine only the solution for the forward part of 
the airfoil and uses the method of characteristics to complete the solu- 
tion for the rear of the airfoil. Further work of Gullstrand is pre- 
sented in references 2k and 25 . 

A list of symbols is contained in the Appendix. 

BASIC EQUATIONS ' 


The basic equations necessary for the discussion of inviscid tran- 
sonic flow consist of a set of partial differential equations relating 
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the velocity components and their gradients at every point, together 
with the auxiliary relation giving the velocity jump through a shock 
wave. For thin airfoils inclined at zero or sma ll angles of attack, 
the differential equations can be simplified by assuming that the shock 
waves are sufficiently weak that the flow is irrotational and isentropic, 
thereby permitting the introduction of a velocity potential <I>. The 
quasi -linear partial differential equation satisfied by $ can he 
expressed in the form 


(a 2 - 4 X 2 ) *3^ + (a 2 - $ y 2 ) Oyy + (a 2 - z 2 ) ® zz - 2® x ® y ® xy - 


2®y^z®yz " ~ ® (l) 

where the subscript notation is used to indicate differentiation and a 
is the local speed of sound given by the relation 

a 2 = ao 2 - 2^1 ($ x 2 + $ y 2 + _ Uq 2) (2) 

In this latter equation U Q and ao are, respectively, the free-stream 
velocity and the speed of sound in the free stream, and y is the ratio 
of specific heats (for air, y = 1.4). 

It is convenient to introduce the perturbation velocity potential cp, 

where 


cp = -U Q x + $> (3) 

If it is assumed that all perturbation velocities and perturbation veloc- 
ity gradients (represented by first and second derivatives, respectively, 
of 9 ) are small and that only the first-order terms in small quantities 
need be retained, equation (l) simplifies to the well-known Prandtl- 
Glauert equation of linear theory 


(1 - *b a ). <P« + + <P„ - ° (4) 

where the free-stream velocity is directed along the positive x axis 
as shown in sketch (a) and where Mq = U 0 /a 0 is the Mach number of the 
free stream. It is well known that equation ( 4 ) leads to useful results 
in the study of subsonic and supersonic flow about thin wings and slender 
bodies but that it is incapable, in general, of treating transonic flows. 
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The failure of linear theory in the 
transonic range is evidenced "by the 
calculated values of <P X growing to 
such magnitude that they can no 
longer he regarded as small quanti- 
ties when compared with U Q . 

Second-order theory for thin 
wings would involve solution of the 
equation 


Sketch (a) 


(See ref. 26, p. lto.) Actually, we are interested in retaining higher- 
order terms only to the extent that is necessary to allow the study of 
transonic flow. Examination of the known properties of transonic flow 
fields indicates that the first term on the right can often become of 
importance and should be retained. It is assumed in the small- 
disturbance theory of transonic flow (refs. 1 through 25), however, 
that the remainder of the terms on the right can be safely disregarded. 
The simplified equation is 

(1 - M/) ¥« + *„ ♦ V zz - Mo 2 ^ = * VW 

where 

k = Mo 2 7T (T) 

u o 

As a result of minor differences in the perturbation analysis, recent 
papers have used at least three other expressions for k. This point 
will be discussed further at the conclusion of the present section. 

Equation (6) is valid only in regions where the necessary deriva- 
tives exist and are continuous. Since these conditions do not hold 
where shock waves occur, an additional equation is needed for the tran- 
sition through the shock. The necessary equation is provided by the 
classical relation for the shock polar (e.g., ref. 26, p. 108) . 


= Mo" 9x*xx + M V + ^z) + 

~ ( 'PyT’xy + ^z^xz^ J 
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'Kt O ~ O 

ir 1 c - 


+ w b 2 = (u a - u b ) 2 


~ ^ „p 

u a u b " a * 


2 a. 2 ~ ^ jc.2 

— u a - u a u b + a* 


(8) 


where u, v, and w refer to Cartesian velocity components with u 
being parallel to the flow direction ahead of the shock, the subscripts 
a and b refer to conditions ahead of and behind the shock, and a* is 
the critical sound velocity, which can be expressed in terms of U 0 
and Mo as follows: 


^ = / 7-1 | 2 

U 0 J 7+1 Mq 2 (7+1) 


(9) 


The appropriate simplified equation is obtained from equation (8) by 
resolving the velocities into components parallel to the axes of the 
coordinate system and carrying out a small -perturbation analysis analo- 
gous to that performed in the derivation of equation (6). In this way, 
the following relation is found between the perturbation velocity com- 
ponents on the two sides of the shock wave: 

(1 - Mq 2 ) (u a - u b ) 2 + (v a - v b ) 2 + (w a - wt,) 2 = Mq 2 (u a -u b ) 2 

Uq \ 2 / 

= k (^7^) ( U a- U b) 2 (10) 

where u, v, and w are the perturbation velocity components parallel 
to the x, y, and z axes. This equation corresponds to the shock-polar 
curve for shock, waves of small strength inclined at any angle between 
that of normal shock waves and that of the Mach lines. On either side 
of the shock wave, the perturbation velocity components are related to 
the perturbation velocity potential in the usual manner 


u 



dcp dcp 

8? w = 81 


( 11 ) 


In addition to satisfying the differential equation and the shock- 
wave equation, the perturbation potential must provide flows compatible 
with the following physical requirements: (a) the perturbation veloci- 
ties must vanish far ahead of the wing and (b) the flow must be tangen- 
tial to the wing surface . Therefore, the following boundary conditions 
are to be specified for the perturbation potential: 
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at x = - w 



at the wing surface, W 


(12) 


_i_ /cwA = ( 13 ) 

u 0 \^z/ w ^x 


where dz/dx is the local slope of the wing surface in the x direc- 
tion. Furthermore, it is consistent with the assumption of small dis- 
turbances to satisfy the second boundary condition on the two sides of 
the xy plane rather than on the actual wing surface. Equation (13) 
thus becomes 


* 

Uo \dz/y 

where the shape of the 


= M = T d x^ y\ 

U 0 \Sz / z=0 dx d(x/c) V c 

wing profile is given by 


£ = Tf 


c 



(no 


(15) 


In addition, it is presumed necessary to prescribe that the direct 
influence of a disturbance in the supersonic region proceeds only in 
the downstream direction and that the Kutta condition applies whenever 
the flow velocity at the trailing edge is subsonic. 

Upon solving the above boundary-value problem for the potential, 
one may determine the pressure coefficient for planar systems by means 
of the formula 



2_ ckp 

u 0 5 x 


(16) 


COMPARISON WITH OTHER STATEMENTS OF THE TRANSONIC -FLOW EQUATIONS 


As a result of minor variations in the perturbation analysis, recent 
papers have used at least four different relations for k, the coefficient 
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of the nonlinear term in the simplified equation for the perturbation 
velocity potential. As indicated in the preceding paragraphs, straight- 
forward development of second-order theory leads to the relation 


k = Mo 2 Zii 
Uo 

This is sometimes simplified (e.g., refs. 8, 10, and 12) to 


(7) 


k 


7+1 

Uo 


(17) 


by arguing that Mq can be set equal to unity in this term without any 
loss in accuracy since the right-hand side of equation (6) is merely an 
approximation to allow the treatment of transonic flows and rapidly 
diminishes in magnitude as Mo departs from unity. In some treatments 
(e.g.,. refs. 7 and 23 ) , equation (l) is divided by a 2 and the quo- 
tient l/a 2 in each term is expanded in a binomial series. When this 
is done, the coefficient k of the term involving cp x cp^ is 


k = Mq 2 


2+ (7-1)Mq 2 ~ 
Uo 


(18) 


Still another expression for k is used by Oswatitsch in the papers 
that form the principal references for the present work. Two deriva- 
tions are given, one based on mass -flow considerations (ref. 21) and 
the other (ref. 20) on simplifying equation (l) under the assumption of 
nearly parallel flow to 


(l - M 2 ) +9yy 

and substituting the following series 


+ cp zz = ° (19) 

for the variable coefficient (l-M 2 ) 


1 - M 2 = 1 - Mq 2 


1-Mq 2 8 cp 
a* - U Q dx 


( 20 ) 


where M is the local Mach number and a* is the critical sound veloc- 
ity as defined in equation (9). Comparison of equations (19) and (20) 
with equation (6) shows that the coefficient k in this approximation 
is 


k = 


1 -M 0 2 
a* - U Q 


( 21 ) 
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A similar situation arises in the derivation of the simplified 
equation for the shock polar. Here again the precise form of the expres- 
sion for the coefficient k of equation (10) depends on the details of 
the perturbation analysis. The most important point from the present 
point of view is that the same expression for k is used in both the 
equation for the potential and that for the shock polar, namely, equa- 
tions (6) and (10). While this point has not always been explicitly 
stated, it is actually a necessary condition for the existence of the 
well-known transonic similarity rules. 

The foregoing discussion has been based on equations obtained by 
assuming that the local velocities are only slightly different from the 
free -stream velocity. On the other hand, many of the recent papers on 
transonic flow about wings and bodies have been based on equations 
obtained by assuming that the local velocities are only slightly dif- 
ferent from the critical sound velocity a*. (See refs. 1 through 6, 9r 
and 13 through l8.) It is shown in reference 12, however, that the 
pressure, force, and moment results obtained using these equations are 
identical with those obtained using the present equations if k is 
selected as given in equation (17)* These results, however, can be 
easily converted to those that would be obtained using any of the other 
expressions for k by simply replacing (7+l)/u 0 by k wherever it 
occurs. 

In order to facilitate comparison with previous results and to 
achieve an economy of notation, the present analysis is carried as far 
as possible without specifying a particular relation for k. That is, 
the equations of the analysis and the reduced parameters with which the 
results are expressed are written containing k which may be equated 
to any of the four stated expressions. However, the actual values of 
the pressure coefficient and Mach number for an airfoil of specific thick- 
ness ratio depend on which relation is selected for k. The present 
calculations have been made using the expression -for k given in equa- 
tion (7)* The principal reason for this choice is that it appears to 
provide a set of equations, or a mathematical model, which approximates 
certain essential features of transonic flow with superior accuracy. 

Before proceeding with discussion of this point, it should be noted 
that the four alternative expressions for k are identical for M 0 = 1, 
and all but that given by equation (17) are zero for Mq = 0. 

A significant case where the four relations lead to different 
results is the prediction of the variation with free-stream Mach number 
of the critical pressure coefficient C D , defined as the value of the 

pressure coefficient Cp at a point where the local Mach number is 
unity. It is important that a reasonably good approximation be main- 
tained for the variation of C ^ with Mo because shock waves make 

F cr 
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their first appearance and the airfoil first experiences a pressure 

drag when C n becomes more negative than C-d somewhere on the air- 
p F cr 

foil surface. In the present approximation, Cr> corresponds to that 

F cr 

value of Cp and, hence, of cp x > at which equation (6) changes locally 
from elliptic to hyperbolic type. This condition is recognized by the 
vanishing of the coefficient of cp^, thus, 

1-M 0 2 - k(9 x ) cr = 0 
or, in view of equation (l6) 


' p cr u Q ^ x ^cr 


(1-Mo 2 ) 


The exact relation for isentropic flow is (e.g., ref. 27, p. 28) 


— + 21li Mq 2 ) 7 -i . 2 . 
7+1 7+1 / 


The variation of 


Cp cr with Mq has "been computed using the exact rela- 


tion and each of the four approximate relations. The results are pre- 

I ‘* 0 1 1 : — sented graphically in sketch (b). It 

\ may be seen that a reasonably good 


i 



Sketch (b) 


approximation for 


is obtained 


over a wide Mach number range when k 
is taken as Mo 2 (7+l)/U 0 or 
(1 - Mq 2 ) /(a*-Uo) , and that a somewhat 
greater error is incurred when k is 
equated to Mq 2 [2 + ( 7 -l)Mo 2 ) /u q . On 
the other hand, a very poor approxima- 
tion results if k is equated to 

(7+1)Aj 0 . 

Similar comparisons can be made 
for local Mach numbers M other than 
unity by noting that the coefficient 
(1 -Mq 2 - kcp x ) of cp^ in equation (6) 
corresponds , in the present approxima- 
tion, to 1 - M 2 , thus. 


1 - W = 1 


kcp x = 1 “ Mq 2 + C, 
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The corresponding exact relation for isentropic flow is 



The results so obtained are generally similar to those indicated in 
sketch (b) although the relative accuracy of the better approximations 
changes somewhat with the situation. All the approximations are exact, 
of course, when Cp = 0. On the other hand, none of the approximations 
are exact, except for isolated cases, when C p is different from zero, 
even though all of the approximations agree among themselves when the 
free -stream Mach number is unity. In order to provide some information 
regarding the errors that are likely to be incurred when C p is not 
very small, sketch (c) has been pre- 
pared illustrating the variation of 
local Mach number with pressure coef- 
ficient for a free-stream Mach number 
of unity. 

A second case where the exact and 
approximate relations can be compared 
is furnished by considering the veloc- 
ity jump through a shock wave. If the 
flow ahead of the shock wave is uni- 
form and parallel to the x axis, the 
results may conveniently be repre- 
sented by the shock -polar diagram in which + 'Wjj 2 is plotted as a 

function of u^. The exact relation is furnished by equation (8). The 
corresponding approximate relations are determined from equation (10) 
by setting u a , v a , and w a to zero, whereby 


v b 2 + w b 2 = [ " ( 1 ” M ° 2 ) + I u b "b 2 ( 26 ) 

Once the variation of (v^ 2 + w b 2 ) with u^ is determined for a given Mq, 
the corresponding variation of (v b 2 + w^ 2 ) with u^ may be readily deter- 
mined since, for this case, 

*b = U o + u b> ^b = v b> *b = w b 



( 27 ) 
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The variation of a/v-^ 2 + w^ 2 with u^ for Mq = 1.2 has been computed 
using both the exact and approximate relations y and the results are 
presented graphically in conventional shock-polar form in sketch (d). 


pf * Sf'o 9 



It is evident from this comparison 
that the best approximation to the 
shock -polar curve is that obtained 
by equating k to Mq 2 (7+1) /u Q . 
Since all shock waves are assumed to 
be normal to the flow direction in 
the course of the present analysis , 
a notable point is that this expres- 
sion for k leads to the exact 
relation for the velocity jump 
through a normal shock wave. 


Sketch (d) 

In problems such as we are considering here, the final test is 
provided by comparison with experimental results. Although both experi- 
mental and theoretical results for the transonic speed range are limited, 
complete information does exist at the present time for the drag of a 
single-wedge section followed by a straight section extending far down- 
stream. (See refs. 13, l4, 1 6, 28, and 29.) The theoretical results 
were determined originally using equations obtained by assuming that 
the local velocities are only slightly different from the critical 
sound velocity a* and are therefore identical with those that would 
be obtained using the present equations, provided k is equated 
to (7+1) /U Q . Sketch (e) shows the theoretical and experimental results 
plotted in the same manner as in the original papers. The small vertical 


(y + l)</3 ~ 

(t/c) 5/3 d o 



Sketch (e) 
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lines on the experimental data points represent the uncertainty of the 
values. This sketch indicates that the theoretical and experimental 
re stilts are only in general qualitative agreement when k = (7+1) /u Q • 

The same results are replotted in sketch (f) with k equated 
to Mo 2 (7+l)/u 0 rather than (7+1) /u 0 . It can "be seen that the theo- 
retical and experimental results are now in nearly perfect agreement. 

(t/c) 5 ' 3 d ° 



Sketch (f) 

Comparison of sketches (e) and (f) provides striking evidence supporting 
the contention that k should he equated to Mo 2 (7+l)/u 0 rather than 

(7+1) /u G . 


DERIVATION OF INTEGRAL EQUATIONS FOR TRANSONIC FLOW 


In order to make the present work more self-contained, a deriva- 
tion of the integral equations for transonic flows having subsonic free- 
stream velocities will be presented even though this has been done 
previously by both Oswatitsch and Gullstrand (refs. 20, 21, and 22). 

The present derivation, in common with that of Gullstrand, proceeds 
through the application of Green' s theorem in a manner closely analogous 
to that employed in linearized wing theory (e.g., ref. 30), except that 
proper cognizance must be taken of the shock discontinuities and of the 
additional nonlinear term in the differential equation for the perturba- 
tion velocity potential. For the sake of completeness and to illustrate 
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the simplifications introduced by making additional restrictions and 
assumptions, the derivation will be carried through for lifting wings 
of finite thickness and span, even though the applications contained 
herein will be confined to two-dimensional flow about symmetrical air- 
foil sections at zero angle of attack. 

The differential equation fundamental to the following discussion 
is equation (6). 


S 2 <p S 2 cp S 2 cp Sep S 2 cp 

(l - M„ 2 ) + + = k 

Sx 2 Sy 2 Sz 2 Sx Sx 2 


(6) 


Since the principal object of the following analysis is to determine 
the pressure which, according to equation (l6), is linearly proportional 
to the perturbation velocity component u, it is convenient to work with 
an equivalent equation for u obtained by differentiating equation (6) 
with respect to x; it is 


(1-Mq 2 ) 


d 2 ^. S 2 u + ^ d 2 / u 2 \ 

Sx 2 Sy 2 Sz 2 8x 2 \ 2 / 


(28) 


It is advantageous to normalize the equations by letting 


x = x, y = 3y, 


z = 3z, 



Sx 



v 


Sep 

dy 



S<£ = k_ 

z p 



(29) 


where 

3 = Vi-Mq 2 

In this way, equations (6) and (28) reduce to the following: 

Sfcp + Sff + - _ 89 S^P 

Sx 2 Sy 2 Sz 2 Sx Sx 2 


(30) 


¥1 + ^ + ^ = v 2 u 

Sx 2 Sy 2 Sz 2 


u 2 


_df_ 

Sx 2 


(31) 
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Before proceeding, it should be noted that the introduction of the 
reduced perturbation velocity component u permits the ready recogni- 
tion of regions of subsonic and supersonic velocities and emphasizes 
the points at which sonic velocity occurs. This relationship becomes 
immediately apparent upon substituting the definition of u into 
equation (24). Thus, 


1 - M g 
1 -Mo 2 


k 

u 

1 - Mq 2 


= 1 


u 


(32) 


from which it is clear, for flows having subsonic free -stream Mach num- 
bers (Mq < l) , that u < 1 when the local velocity is subsonic, u = 1 
when it is sonic, and u > 1 when it is supersonic. 

As noted above, the derivation proceeds from Green's theorem which 
relates a volume integral over a region R to a surface integral over 
the surface £ enclosing R. If o and Q are any two functions which, 
together with their first and second derivatives, are finite and single- 
valued throughout R, Green's theorem states 


//(" £ - 0 l) 4£ ’ -fff ^ “ <33) 

E R 

where the directional derivatives on the left side are taken along the 
normal n, drawn inward, to the surface E. It is convenient to let 
G = u and to choose 0 as the fundamental solution l/r 3 of the 
equation V 2 o =0 



1 


(x - l) 2 + (y - ti) 2 + (z - l) 2 



whereby equation ( 33) becomes 


(34) 


_L M _ -a A ( ’JLY 

r 3 <3n dn \ r 3 /_ 

R R 

The variables of integration in the equation are I , fj, C while x, y, 
are the coordinates of a point P. It must be observed that l/r 3 is 
singular at r 3 = 0 and u is discontinuous at the shock wave. The 
point P and the shock wave must, therefore, be excluded from the 
region R. The exclusion of P from the region R is accomplished by 
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Sketch (g) 


enclosing it within a small sphere. 

The shock wave is excluded by altering 
the boundary of the region so that it 
goes around the shock wave. In this 
way, equation (35) may be applied to 
the region R u bounded by the xy 
plane and a hemispherical dome of infi- 
nite radius lying above this plane, 
exclusive of the subregions surrounding 
P and the shock wave (see sketch (g)). 
Since, furthermore, the values of u 
may be assumed to diminish sufficiently 
rapidly with distance that the contri- 
butions of the integrals over the hemi- 
sphere vanish, the following result is 
obtained: 


u(x,.y,£) = - -L j 

4* J 


+ 00 \ - 

1 ^ u ' 


du u _ a/iA] 
hi UU ht\r 3 )_ 


dfdfj 
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— (^)dR 
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(36) 


where the subscript u denotes conditions on the upper side of the xy 

v the subscri P ts a and denote values immediately ahead of and 

behind the shock wave, and S is the surface of the shock wave. The 
volume integral is defined as follows when P is ahead of S. (For 

sake of brevity, \|r is written in place of — -§f- ( \ 

r 3 M 2 \zj 
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(37) 
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It is clear that the corresponding definition of the volume integral 
when P is behind S can be obtained by rearranging the limits of 
integration. If P is kept fixed in the upper half space and the 
region R j bounded by the xy plane and a hemispherical dome of infi- 
nite radius lying below this plane is considered, it follows in a simi- 
lar manner that 


0 = 4 Mk ir - (^)] a!d * ' 4//{fc £ ■ - 6 i (4)] 

Sj. £ 


it 


— _ u JL (jL 


dn 


— ( - — ) | ds - — rrr ~ — 

Bn ' r 3 < _J^J kit JJJ r 3 \ 2 


-2 \ 

)d R 


(38) 


where the subscript 2 denotes conditions on the lower side of the xy 
plane and the volume integral is defined as follows: 

Mkw$)'* a lIF* R ‘&o /”*/*{/* id|+ / wC } 

R z R l %+ e 2 

(39) 


Introducing the notation 


Au = u u - u 


± Bu ^2 

l ’ B£ BC BC 

and adding equations ( 36 ) and ( 38 ), we have 
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where the surface integral over S represents the sum of the corre- 
sponding integrals over S u and Sj and the volume integral over R 
represents the sum of the integrals over R u and R}. The integrand of 
the first integral of equation (^l) is zero over all of the xy plane 
except the surface W of the wing sufficiently extended to include the 
effect of the edge singularities and is, in certain cases, exactly equal 
to the value of u given hy the linearized theory of subsonic flow 
about thin wings (e.g., ref. 30 ) * 



W 


It can be seen that the first integral of equation (4l) may be equated 
to u^ when the problem is one in which 

Au = Au l 

and 

A(5u/^C) = A(Su/^C) l 

This condition exists in those problems where Au and A(c)u/ 5£) are pre- 
scribed at the outset by the boundary conditions; for example, (a) given 
the loading on a lifting surface, find the camber distribution; 

(b) given the shape of a symmetrical nonlifting airfoil, find the pres- 
sure distribution. 

Equation (kl) can be regarded as the final integral equation for u, 
but it is advantageous for the forthcoming analysis to perform two more 
operations. They are to integrate the volume integral twice by parts 
with respect to I, taking proper cognizance of the definitions given' in 
equations (37) and (39)> and to decompose the surface integral over the 
shock wave into components parallel to the axes of the coordinate system. 
In this way, the following equation is obtained: 
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Although the integration by parts of the triple integral performed 
in going from equation (4l) to (43) may seem somewhat arbitrary, the 
resulting equation is superior from the point of view of obtaining 
approximate solutions. For example, the triple integral of equation (4l) 
shows a very strong influence of the velocities in the region immediately 
surrounding P since they are multiplied by l/r 3 . This influence is 
largely nullified in the triple integral of equation (4-3) because part 
of the region has a negative influence and part has a positive influence. 
The predominent influence in the latter case is furnished by the term 
u 2 /2 standing outside the integral. The contribution of distant regions 
is also diminished in importance in the triple integral of equation (43) 
since their influence varies inversely with the third power of the dis- 
tance, rather than the first power as in equation (4l). A further 
advantage is that the value of the triple integral of equation (43) is 
continuous through a shock wave rather than discontinuous as is the case 
with equation (4l). A point of great importance in the approximate 
solution described herein arises from the fact that the integration by 
parts provides extra terms (those containing u 2 /2) in the integrals 
along the shock surface S which combine with those already present in 
such a way that the contribution of these integrals becomes very small 
when the shock waves approach normal waves, as is usually the case at 
high subsonic speeds. 
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In addition to satisfying the integral equation for u given in 
equation (43), the velocity components on opposite sides of shock waves 
must be in accord with the simplified relation for the shock polar given 
in equation (10). This equation may' be rewritten in normalized form by 
introducing the quantities defined in equation ( 29 ), thus 


(u a - u b ) 2 + (v a -v b ) 2 + (w a -w,) 2 = 


(u a - u b ) : 


'(W) 


Two alternative forms of equation (44) are the following: 

(y 1 ' q ) ^ u a -u b) 2 + ( v a " v b) 2 + ( w a - "b') 2 = 0 (^5) 

and 


(u a - u b ) 



+ ^a - ^b) 2i + '(w a - w b ) 2 = 0 


(46) 


If the shock wave is a normal wave and the flow. is parallel to the x 
axis (i.e., v a " = ^a ~ ^b = "but u a 4 > it can be seen from equa- 

tion (45) that the normalized perturbation velocity component u jumps 
from 1 + A immediately ahead of the shock to- 1 - A immediately 
behind , the- shock • On the other hand, equation 1 . (46) shows that the quan- 
tity u - u 2 /2 is, equal on the two sides of the shock. This is con- 
sistent with -the fact that the latter quantity corresponds, in the 
transonic approximation, to the mass flow, which is continuous through 
a normal shock. - 

’•The solution, by the present method, of the general problem of 
transonic flow about thin wings requires the solution of equation (43) 
while taking proper, account of the shock relations .given in equation (44). 
This .represents a formidable task well beyond the reach of the present 
analysis.- Simplification can be achieved in two. ways: by restricting 

attention to a less general class of problems and by introducing addi- 
tional simplifying approximations. The first way is, of course', much . 
to be preferred. Accordingly, in most of the following analysis, atten- 
tion will be confined to two-dimensional flows. The necessary equations 
can be obtained from equations (42) through (44) above by integrating in 
the f\ direction from fj •= - 00 to fj = +», noting that v = 0* and that 
u and w are independent of fj. They *are as follows: 
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u 



In — Au In 

r 2 ac 


* f - tM ^ i 



where 

_1_ 1 

r 2 f . 2 -|l/2 (^9) 

(x - f ) 2 + (z - l) J 

and 


(u a -u b ) 2 + (w a -w b ) 2 = ( — g Ub ) (u a -U b ) 2 (50) 

As remarked following equation (42), considerable simplification 
results in both two- and three-dimensional problems if attention is 
confined to the determination of the pressure distribution on symmetri- 
cal nonlifting wings of specified geometry. This restriction permits 
the introduction of the relations 


Au = Au l = 0 

A | = A |* A (|) L a8lven (51) 

into the integral over W in either equation (43) or (48). This inte- 
grai is then equal to the linear-theory solution u L given in equa- 
tion (42) or (47) and can be determined completely at the outset of the 
analysis. 
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SIMPLIFICATION AND APPROXIMATE SOLUTION OF INTEGRAL EQUATIONS 
FOR TRANSONIC FLOW ABOUT NONLIFTING WINGS 

Qualitative Discussion of Integral Equations 


The integral equations and the auxiliary relations developed in the 
preceding section provide a means for the discussion of the aerodynamics 
of symmetrical nonlifting thin wings of specified geometry in flows with 
free-stream Mach numbers up to unity. At the present stage of develop- 
ment, however, it is necessary to introduce some further approximations 
before solutions can be obtained. One of the more accurate of these 
involves two statements about the nature of the shock waves. They are: 

(a) all shock waves are assumed to lie in a plane perpendicular to the x 
axis, and (b) the shock waves are assumed to be normal shock waves (i.e., 
normal to the local flow direction). These two statements are slightly 
contradictory in themselves but might be expected to approach the true 
conditions quite closely for flows about thin airfoils. The first state- 
ment corresponds to setting cos(n,f}) and cos(n,C) to zero, thereby 
eliminating part of the integrals over S of equations (43) and (48). 

The second permits an advantageous introduction of equations (30) an< 3- (46) 
to eliminate the remainder of the integrals over S. The above assump- 
tions correspond to setting v and w to zero both before and after the 
shock wave and lead to the following relations: 



Equations (43) and (48) thereby simplify to 



R 


for three-dimensional flows, and 



for two-dimensional flows. These equations correspond to those used by 
Oswatitsch and Gullstrand (refs. 20, 21, and 22) although a number of 
further assumptions were necessary before approximate solutions could be 
obtained for the velocity distribution on an airfoil surface. The pres- 
ent analysis also requires many of the same or similar assumptions, but 
there are a number of general points which should be discussed before 
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any further approximations are introduced. These points were not 
observed in the previous papers and the results suffer by lacking unique- 
ness in some cases and failing to converge in other cases. 

Before proceeding further, it should be observed that the solutions 
of equations (53) and (5*0 must approach those of linear theory when the 
free-stream Mach number is much less than unity, since u « 1 and the 
terms involving the square of u become negligible with respect to 
those linear in u, thereby leaving only 

(u) ^ u T 

v 'Mo« 1 L 

In the interest of brevity, the integral equations for u given in 
equations (53) and (5*0 may be rewritten as follows: 



where 



(56) 


1 


r_L r 


LwJ 

7 2 


R 


for three-dimensional flows, and 


dldfjdC 


(57) 


I 



(58) 


for two-dimensional flows. Although I is a function of u and is 
therefore unknown, it is informative to rewrite equation (56) by solving 
for u in terms of I and u^, thus 


u = 1 ± </ I - (2u L -l) = 1 ± a/ I - L 


(59) 


where 

L = 2u L - 1 

Several points are to be observed at once with regard to equation (59)* 
First of all, the discriminant must always be positive in order to 
obtain real values for. u, thus 


I > L 


(60) 
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Furthermore, the choice of the plus or minus sign determines whether 
the local velocities are subsonic or supersonic. A change in sign at 
a point where the radical is zero corresponds to a smooth transition 
through sonic velocity. A change in sign at a point where the radical 
is not zero corresponds to a discontinuous jump in velocity. As pointed 
out following equation (46), such discontinuities correspond to normal 
shock waves and are permissible when they proceed from supersonic to 
subsonic velocities (or from plus to minus sign in equation (59)) when 
progressing in the flow direction. Discontinuities in the reverse 
direction are inadmissible since they correspond to expansion shocks, 
an inqpossible phenomenon which violates the second law of thermodynamics. 

The values of u^, and hence L, can be calculated for any given 
wing and are generally characterized by certain regions in which u^ 
is positive and other regions in which it is negative. The absolute 
values increase continuously with increasing Mach number and the max imum 
positive values may considerably exceed unity as sonic velocity is 
approached in the free stream. Not very much can be stated at this 
point about the values of I, except that they depend on the distribu- 
tion as well as magnitude of u and that the above inequality must be 
satisfied. The relation between the two curves is of utmost importance, 
however, and will be discussed qualitatively in the following paragraphs. 

In order to remove unnecessary complications and to facilitate the 
discussion, the following remarks will be confined to the relations 
between the functions I, L, u, and u L evaluated at the airfoil surface 
at a single spanwise station situated in the plane of symmetry of the 
wing. In this way, each of the four functions reduces to a function of 
a single variable x and can be illustrated simply by curves rather 
than surfaces or hypersurfaces . The subscript W is appended to 
u and u L to denote that the values are those at the wing surface. In 
that which follows, the curves will be shown on two separate plots, one 
containing the L and I curves representing the components involved in 
the solution of equation (59) and the other containing the u^ and u,- 
curves representing, respectively, the velocity distributions given by 
transonic theory and by linear theory. In order to make the discussion 
more definite, the curves will be drawn qualitatively as they would 
appear for a circular-arc airfoil having its maximum thickness at the 
midchord position. A quantitative discussion of these characteristics 
will be taken up for the same airfoil in a later section. 

The linear-theory solution for subsonic flow about circular- 

arc airfoils can be easily derived through application of the expression 
given in equation (47). It is found that the values of u L are symmet- 
rical about the midchord position at all free -stream Mach numbers less 
than unity. It follows directly that the L curves possess the same 
symmetry. 
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For pure subsonic flows about such an airfoil, it is well known 
that the more exact treatments, such as the Karman-Tsien or the classi- 
cal iteration methods reviewed in reference 19 , indicate that the 
nonlinear-theory solutions and, hence, the u^ and I curves, are also 
symmetrical about the midchord position. Since sonic velocity is not 
attained at any point, the L and I 
the curves for this condition are 
shown in sketch (h) . The second 
part of this sketch illustrates the 
fact that u w possesses larger 
values over the middle of the air- 
foil chord than does uj^. This 
relation has not been deduced from 
the foregoing preliminary consid- 
erations but is known from many 
sources in classical subsonic 
theory. The same result has been found once again in the present work 
by carrying out an approximate solution of equation (59). The details 
of these calculations will be described in a later section. 


curves never touch. Sketches of 




Sketch (h) 


The curves shown in sketch (h) are typical of those for all Mach 
numbers less than the critical Mach number Me r , defined as the lowest 
free -stream Mach number at which sonic velocity (u^ = l) occurs some- 
where on the wing surface. The aforementioned symmetry properties are 
preserved for all Mach numbers less than the critical, but the amplitude 
of all four curves increase with increasing Mach number. For the symmet 

rical circular -arc airfoils considered here, the maximum values of u T 

• W 

and Uy and, hence, L and I occur at the 50-percent-chord position 
for all subcritical Mach numbers. Thus, in addition to the . requirement 
that I > L at every point, it is necessary when the Mach number is 
equal to the critical value, that I = L at the 50 -percent -chord posi- 
tion. A rather interesting additional requirement that follows from the 
quadratic nature of equation ( 56 ) together with the assumption that u 
reaches a smooth maximum at the 50-percent-chord station is that the 
I and L curves have not only the same first derivatives at this station 
but also identical second derivatives. A typical set of curves for this 
Mach number is shown in sketch (i). 

It is to be remarked that the minus 
sign is to be used in equation (59) 
for all free -stream Mach numbers 
equal to or less than M cr . 

It is interesting to contem- 
plate the various possibilities 

that may occur when the Mach number = M cr 

is increased beyond the critical. 

Accordingly, let us first consider 




Sketch (i) 
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the implications of assuming that the velocity distribution uy remains 
symmetrical and continuous across the chord and simply increases in 
magnitude with increasing Mach number. (This is, in fact, the only type 
of mixed flov field that the classical iteration methods have indicated, 
but many doubts have been expressed from time to time regarding the 
convergence of the methods when Mq > M cr .) With such an assumed sym- 
metry of the flow, the four curves will appear qualitatively as shown 
in sketch (j). As may be seen, the curves are all rather similar to 

those previously discussed for lower 
Mach numbers. The outstanding dif- 
ference is the relation between I 
and L. When Mq = M cr , the two 
curves are tangent at the 50-percent- 
chord station, and the radius of 
curvature of the I curve at the 
same station is equal to that of 
the L curve. When Mq > Mc r , 
tangency occurs at two points, 
equally spaced before and after the 
50 -percent -chord station, and the sign in equation (59) is to be changed 
to plus over the portion of the chord lying between the two points of 
tangency. In order for uy to attain its maximum value at midchord, as 
shown in sketch (j), the radius of curvature of the I curve must be 
less than that of the L curve at the 50-percent-chord station. At the 
tangent or sonic points, it follows from equation (59) and- the assumption 
of smooth acceleration or deceleration through sonic velocity that the 
second derivative of the I curve is greater than that of the L curve. 
Similar considerations apply for airfoils that are not symmetrical about 
the midchord station. The occurrence of shock-free flow would again 
require that the I curve be tangent to the L curve at two points along 
the chord. 




M 0 > M cr , Symmetrica! Solution 
Sketch (j) 


Subsonic 


Before leaving this subject, it is interesting to inquire what the 
result would be of a slight alteration of the airfoil shape, assuming 

that the original shape was such 
/' x that the associated flow was of the 

shock-free mixed type. Consider, 
for sake of simplicity, that the 
basic airfoil is symmetrical about 
the 50-percent-chord station. The 
assumed shock-free flow is there- 
fore symmetrical about the midchord 
station and appears as shown in 
sketch (k). In order to preserve 
the geometrical symmetry, consider 
that the airfoil shape is changed by 
the addition to both upper and lower 


U 0 



Sketch (k) 
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surfaces of small lumps located at the 50-percent-chord station. Since 
the flow adjacent to the humps is supersonic, their disturbance pattern 
is propagated downstream in narrow bands which reflect alternately from 
the sonic line and the airfoil surface as shown in sketch (Z). Busemann 


has pointed out in reference 31 that 
of the supersonic region where the 
local Mach numbers of the original 
flow approach unity, the reflections 
become increasingly concentrated, 
the disturbances amplify, and 
finally a shock occurs. It thus 
appears that there are only a 
restricted number of shock-free 
mixed flows and that, in general, 
the supersonic region terminates 
with a shock wave. 


as the disturbances reach the rear 



On the other hand, if the dis- 
turbances were propagated forward they would culminate in a shock wave 
at the forward sonic point. This shock, however, would be an expansion 
shock in which the velocity jumps from subsonic to supersonic values. 

As mentioned previously, expansion shocks are forbidden by thermodynamic 
considerations and, hence, must be excluded from the present analysis. 
This can be accomplished by stipulating that the transition from subsonic 
to supersonic velocities at the forward sonic point be smooth, or that 
the L and I curves retain one point of tangency. In addition to pre- 
venting the occurrence of expansion shocks, it appears, on the basis of 
the foregoing, that this requirement also effectively rules out all the 
undesired forward propagating disturbances. 


The preceding discussion provides an insight into the mechanism for 
the development of asymmetrical flows about symmetrical airfoils and for 
the occurrence of shock waves. In the present work, these considerations 
are reflected in the relation between the I and the L curves. Thus, 
recall that if the flow is shock free as shown in sketch (k) , the I and 
the L curves are both symmetrical as shown previously in sketch (j)# 

If the airfoil shape is changed in the manner indicated in sketch (z), 
it is evident that both the I and the L curves will become altered. 

In so doing, the L curve remains symmetrical about the midchord sta- 
tion, but the asymmetrical nature of the u distribution causes the 
I curve to take on larger values over the. rear of the airfoil than over 
the front. If the flow adjusts itself so that the I curve is tangent 
to the L curve at a point on the forward half of the airfoil, as it 
must do to avoid the occurrence of forbidden expansion shock waves, it 
is likely that the curves will not be tangent at a second point along 
the rear half of the airfoil. In the application of equation ( 59 ) > 
therefore, the sign changes from minus to plus at the point of tangency 
but may change back to minus at a point where the two curves have 
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different values. As noted previously, such a condition corresponds, 
in the present analysis, to the occurrence of a shock wave and is asso- 
ciated with a u w curve of the type shown in sketch (m) . 




M.>M er , Asymmetrical Solution 


Since it has been indicated that 
the shock-free mixed flow is the 
exception rather than the rule, the 
curves of sketch (m) may be regarded 
as typical for the supercritical Mach 
number range. In any event, these 
curves may be considered as the more 
general ones since they include the 
symmetrical curves as a special case. 
These matters will arise again and be 
the topic for further discussion in 


Sketch (m) 

the next section in which an approximate method for the solution of the 
integral equation for transonic flow will be described. 


Simplification of the Integral Equation 


The remainder of the present discussion will be concerned with two- 
dimensional flow about nonlifting symmetrical airfoils of specified 
geometry under the assumption that any shock waves which may be present 
are normal shocks situated perpendicular to the x axis. The analysis 
will be based therefore on equation (5*0 which, when written in full, is 


tshh «x-C) 2 -(z 4) 2 ] dCd £ 

2 [(x-C) 2 +(z-£) 2 ] 2 

(61) 

Approximate solutions of this equation could conceivably be worked out 
numerically by_starting with a two-dimensional grid of suitably selected 
values for u(C, £) and iterating until convergence is obtained. Such 
calculations would proceed by inserting the assumed values for u into 
the double integral and solving to obtain the next approximation for 
u(x,z), making use of the tangency condition on the surfaces or functions 
represented by I and L as discussed in the preceding section. If the 
first approximation for u is taken to be the results given by incom- 
pressible or by linearized compressible flow theory, as in the Rayleigh- 
Janzen and other classical iteration methods, it seems to be the prevail- 
ing belief that convergence will be obtained only when the free -stream 
Mach number is sufficiently small that the flow is subsonic at every 
point. It is at this point that Oswatitsch (refs. 20 and 2l) supplied 
the important idea that mixed flow fields containing shock waves can be 
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obtained if the starting u distribution contains shock waves. Thus, 
in place of taking the starting solution to be that given by incompres- 
sible or linearized compressible flow theory, the idea is to start with 
a reasonable guess for the values of u, being sure to include a proper 
discontinuity complying with the shock relations of equation (52), and 
then to proceed to the solution. As will be seen in the succeeding 
sections, it is not necessary to be highly accurate in the initial 
guess for u. 

A source of difficulty in the numerical solution of equation (6l) 
by an iteration process is the double integral. If it could be reduced 
to a single integral by introducing a suitable approximation, the entire 
problem would be greatly simplified. In the present analysis, it is 
assumed, following Oswatitsch, that approximate knowledge of the velocity 
distribution is sufficient for providing a working approximation for the 
double integral. In particular, it is assumed that a sufficiently good 
approximation to the velocities in the vicinity of the wing can be 
expressed in terms of the local coordinate z, the ordinates of the 
airfoil surface z(x) , and the desired but unknown velocity distribu- 
tion uy(x) on the airfoil surface. This permits one integration to 
be performed, thereby reducing the double integral of equation (6l) to 
a single integral. 

A number of statements regarding the variation of u with z over 
the middle portion of the airfoil can be made immediately. For example, 
u starts from the value u w at the airfoil surface with an initial 
rate of change given by the irrotationality condition 


( 8 ), -( 8 ), 

and probably vanishes at great distances as l/z 2 . These conditions, 
of course, are not sufficient to determine completely the variation of 
u with z, but may be used as the basis for the development of an approxi- 
mate relation. Oswatitsch (refs. 20 and 21) has already considered this 
step and has suggested the following relation: 


u(x,z) 


[l + (z/b)] 2 


(63) 


where b is a function of x so chosen that the irrotationality con- 
dition is fulfilled at z = 0. Thus, differentiation of equation (63) 
with respect to z and insertion of the definitions of equation (29) 
and the boundary condition of equation (l4) yield the following: 



30 


NACA TN 3096 


2 % 

(du/^z) w 


gp 3 u w 

(dw/bx) w k (^w/^x) w 


2p 3 uw _ _ 2uw 

kUo (d 2 Z/<3x 2 ) (S 2 z/Sx 2 ) ' 


(64) 


where Z represents the reduced ordinates of the airfoil, related to 
the actual ordinates by 

kU 

Z(x) = Z(x) (65) 

3 

Attention is called to the fact that the approximate relation for u(x,z) 
given above is not entirely satisfactory. Evidence of this is provided 
by the fact that u is indicated to be zero in the region ahead of the 
leading edge and behind the trailing edge where b is infinite and that 
the discontinuities in u at the shock surfaces are consistent with the 
shock relations only at the surface of the airfoil. The errors in the 
pressures on the airfoil surface resulting from the former are small due 
to the attenuating influence of distance, and those resulting from the 
latter have been partially compensated for by a readjustment of the 
approximation at Mach numbers near unity where the shock strength becomes 
large . 

Substitution of equation ( 63 ) into the double integral of equa- 
tion ( 6 l) permits integration with respect to Thus, by performing 

this integration and setting z = 0 , the following approximate integral 
equation is obtained for u^: 

¥-/¥<¥)« <«> 

o 

The function E is 

E ( E(X) = - | X | ( 5-10X 2 + X 4 ) - 

V V *(i+x 2 ) 5 l - 2 

(1-10X 2 + 5X 4 ) In | X | - — (1+ X 2 ) (25-71X 2 -X 4 - X 6 ) 1 


( 67 ) 
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The nature of E(X) is illustrated 
graphically in sketch (n). Although 
the integration interval is indicated 
in equation (6l) to extend from 
| = -00 to l = +«, the contribution 
of the regions ahead of the leading 
edge and behind the trailing edge is 
zero since b is infinite. The 
integration need, therefore, be 
carried out only over the chord c. 

It should be noted that the inte- 
gral in equation (66) corresponds 
to l/2 in equation (56). 

Although equation (66) is con- 
siderably simpler than equation (6l) 
owing to the replacement of the 
double integral by a single inte- 
gral, many of the essential diffi- 
culties remain since the integral 
equation is still nonlinear and 
the kernel is singular. Since no 
known methods exist for the solu- 
tion of such equations, we can only Sketch (n) 

proceed at the present time by introducing additional simplifications. 
One method proceeds by approximating Uy with some simple functions 
having certain parameters temporarily unspecified. Values for the 
latter are determined by substituting the functions into the integral 
equation and satisfying the equation at a limited number of points equal 
to the ‘number of unspecified parameters. At this point in the analysis, 
Oswatitsch assumed that the variation of u w across the chord could be 
represented by a parabola, one or two half -parabolas, or a rectangle 
combined with a half -parabola, as illustrated in sketch (o), all of 
unspecified height and chordwise 
extent. No account was taken in the 
integral of the influence of the 
region between the leading edge and 
the station where u T = 0 on the 




% 

forward part of the airfoil nor 

between the station where u T =0 

L W 


Sketch (o) 

on the rear of the airfoil and the 


trailing edge. Upon inserting a selected combination of the above- 
mentioned elements into equation (66) and integrating, there resulted a 
system of simultaneous quadratic algebraic equations having as many mem- 
bers as there were elements In the selected general form of solution. In 
many cases, Oswatitsch used only one element, either a parabola or a half- 
parabola and assumed a mean value for b for the entire chord. The 
method included no provision for the improvement of the result through 
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iteration or other means, the only measure of the accuracy being the 
degree of correspondence between the initial and the final uy distri- 
butions. Nevertheless, the results presented in references 20 and 21 
were encouraging in that they showed many of the known characteristics 
of transonic flow about airfoils. In particular, shock waves appeared 
when the critical Mach number was exceeded and moved rearward with 
further increases in Mach number. However, because the initial values 
for Uy were generally substituted not only into the integral but also 
into the term uy 2 /2 standing outside the integral, and because the 
tangency requirement on the I and L curves necessary for smooth tran- 
sition from subsonic to supersonic velocities was not realized, a multi- 
plicity of solutions was often obtained. In one case, three solutions 
were actually given and more could have been obtained which would have 
fulfilled equally well the condition of correspondence between initial 
and final result. In addition, the quality of the results appeared to 
deteriorate when the Mach number was increased to higher values, the ' 
upper limit of acceptability appearing to be a Mach number of about 0.88 
for 6-percent-thick airfoils. 

The integral equation method has been developed further in refer- 
ences 22 through 25 by Gullstrand. In the first of these, reference 22, 
equation (6l) is simplified to a single integral equation through the 
use of a more elaborate velocity-distribution function than that given 
in equation ( 63) , and the resulting equation is solved by an iteration 
process. The introduction of an iteration procedure makes a marked 
improvement over the method of Oswatitsch since it then becomes practi- 
cal to increase greatly the number of elements with which uy is repre- 
sented. The method is applied to three 6-percent-thick NACA low-drag 
airfoils and the resulting velocity distributions are given. In common 
with the original method of Oswatitsch, difficulties occur when the Mach 
number is too close to unity. The highest Mach number for which results 
are given is 0.91* More recently, Gullstrand has presented approximate 
solutions for the velocity distributions on symmetrical airfoils in 
sonic flow in reference 23. 1 In both of these works, however, the itera- 
tion process proceeds by inserting the known values into both the inte- 
gral and the term Uy 2 /2 standing outside the integral. This procedure 
is equivalent to replacing the second-degree equation for uy with a 
linear equation and obscures or loses many of the characteristics of the 
quadratic solution discussed in the preceding section. 


1 Since equation (6l) degenerates at a Mach number of unity where 3=0, 
the sonic results of reference 23 are obtained by first developing a 
new double -integral equation to replace equation (6l) and then simpli- 
fying and solving by an iteration procedure. The integral equation is 
only used, however, to determine the velocity distribution over the 
portion of the airfoil forward of the station of maximum thickness. 

The remainder of the solution Is obtained by means of the theory of 
characteristics . 
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The present analysis also proceeds through the use of an iteration 
scheme, partly numerical and partly graphical, but the known values are 
substituted only into the integral at each step of the process. The 
quadratic nature of the equation is thereby retained and the tangency 
condition on the I and L curves can be fulfilled. Satisfaction of 
the latter requirement is essential for uniqueness and convergence of 
the solution. Inasmuch as the calculations for the circular-arc section 
were well advanced when Gullstrand's papers were received and were pro- 
ducing reasonable results with the simple velocity-distribution function 
of Oswatitsch, it was decided to continue rather than to start over 
using a more elaborate velocity-distribution function. It has been 
found, however, that additional attention must be paid to the influence 
of the region behind the trailing edge at the higher Mach numbers. 

Upon observance of these additional refinements, results are obtained 
for all Mach numbers up to unity. At the lower Mach numbers, these 
results are in general agreement with those found by the simpler, 
although more approximate, method of Oswatitsch. The present method 
carries right on, however, into the higher Mach number range where the 
simpler method met with difficulties and succeeds in showing the well- 
known invariance of local Mach number with changes in the free-stream 
Mach number as the latter approaches unity. 

In the present calculations for circular-arc airfoils, the neces- 
sity for additional refinement begins at Mach numbers somewhat greater 
than that at which the shock wave reaches the trailing edge. Conse- 
quently, no attempt is made to account for the influence of the region 
behind the wing until the Mach number becomes sufficiently large for 
the shock wave to reach the trailing edge. At higher Mach numbers, the 
influence of this region is approximated in the following manner. First, 
it is assumed that the shock wave which stands at the trailing edge 
remains of the strong family as the free-stream Mach number is increased 
to unity. It is assumed furthermore that the flow is parallel to the 
x axis at the shock position and that the shock wave is normal to the 
local flow. With these assumptions, it follows from equation (V?) that 
the values of u immediately behind the shock wave are related to 
those immediately ahead of the shock by 

= 2 - % ( 68 ) 

Since u a is given by equation (63), it follows that u^ can be 
expressed in terms of the values of uy and b immediately ahead of the 
shock wave in the following manner: 


u 


u-u(z) = 2 - 


■w n 


[ 1 + (z/b a )] 2 


( 69 ) 
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Because the shock terminates within a finite distance of the airfoil at 
all subsonic Mach numbers, equation ( 69 ) is appropriate only within a 
certain range of \z\, namely, that for which u a > 1. In the present 
calculations, the end of the shock occurs at such large |z| when the 
shock is at the trailing edge that equation ( 69 ) has been used to repre- 
sent the values for u behind the trailing edge for all |z|. It is 
further assumed that the contribution of the region behind the trailing 
edge can be satisfactorily approximated by equating u to u^ for all 
points behind the trailing edge. It is recognized that neither of these 
assumptions constitutes a good approximation for u at great distances 
from the wing, but the attenuating influence of distance diminishes the 
error in the values of the integral at points on the airfoil surface. 

In this way, equation (66) comes to be replaced with the following rela- 
tion determined from equation (6l) by performing an integration with 
respect to £ under the assumption that the variation of u with z 
(or £) is given by equation ( 63 ) for stations ahead of the trailing 
edge and by equation ( 69 ) for stations behind the trailing edge. 


u w 2 

U W = u Lw + — 




- 2 +oo 

a / _fi~x 


d! + 


Ui 

2b, 


+00 


- f f (— 

ba J c \ ba / 


(70) 


where E is as defined by equation ( 67 ) and F is given by 

£1 


F(Z) 


-/ 



F l ^ ) = F(X) 


it(l+X 2 ) 3 L 


*|x|(3-x 2 ) - 


1.0 1.5 2.0 2.5 3.0 

X 


2(l-3X 2 )ln|x| - ( 3+2X 2 - X 4 ) 

(71) 

The nature of F(X) is illustrated 
graphically in sketch (p) . In this 
case it is apparent that the three 
integrals of equation ( 70 ) taken 
together correspond to l/2 in equa- 
tion (56). 


Sketch (p) 


.5 
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To summarize, equation (66) is used in the present calculations 
when the Mach number is less than that at which the shock wave first 
reaches the trailing edge, and equation (70) is used for higher Mach 
numbers. No significant discontinuities are produced in the value of 
the integrals, however, since the contribution of the additional inte- 
grals only becomes significant at Mach numbers greater than that at 
which the procedure is changed. 


Numerical Evaluation of Integral 


One of the principal steps in 
the iteration method used herein for 
the solution of equations (66) and 
(70) is the evaluation of the inte- 
grals. Since u^ and b are gener- 
ally prescribed by a set of numerical - 
values rather than by analytical * 

functions, a numerical technique has 
been used for the integration. This 
process consists of replacing the 
prescribed u^ distribution with a 
stepwise approximation as indicated 
in sketch (q), introducing a mean 
value for b for each of the rec- 
tangular elements, integrating to 
determine the contribution of a 
single element, and summing the 
influence of all the elements. The contribution of a single element of 
width 2 situated on the airfoil chord, as typified by the shaded area 
of sketch (q), is given by 




X/t 


2 


( 72 ) 
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Upon performance of the indicated operations, the following expression 
is obtained for f 1 : 


A » 


+ '^2)4 jjf TaJ ( 1 + A 2 ) 4 - ( 1+ A 2 ) 2 + 8(1 + A 2 ) -8 +12A(A 2 -l)ln | A | - 


12(1 + A 2 ) 


A(l+ A 2 ) [ (1+ A 2 ) 2 + 12] j + - 

J 12(1 + B 2 ) t 2 


12(1 + B 2 ) 

8(1+ B 2 ) “8 J + 12B(B 2 - l)ln|B| - B(l+B 2 ) 


J3_ 

B 


(1 + b 2 ) 4 -(i+b 2 ) 2 + 


(1 + B 2 ) 2 + 12 


( 73 ) 


where 


A 


= l i +g (x- ^j) = l/2lA l^-IjN /22j\ 

2b t 4 \ hi/ + 2 \ l ± J V hi ) 


B = ~ ~ ^1^ _ 1 f ^ ^i^ j 1 f f ^i ’N 

2bi 4 \bi/ 2\ Zi / V bi/ 


Thus, the integrals in equations (66) and (70) that contain E are 
approximated as follows: 


(74) 


f %•&■)•< -i =M( 


x-ii\ / 2Z^ 




( 75 ) 



It is evident from equation (70) 
that the contribution of a single 
element of width Z situated behind 
a shock at the trailing edge, as typi- 
fied by the shaded area of sketch (r), 
is composed of two parts. The first 
depends on an integral involving 
Tt/c E(x) and is evaluated using fi in 
the manner just described for elements 
on the airfoil. 


Sketch (r) 
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The second depends on an integral involving F(X) 



which, upon evaluation, yields 


l Sl = (^ 2)2 {2 TIT C(l + A2)2 + d + A 2 )-2]-2A In |a| -A(l+A 2 )j + 
— If t|t [(l + B 2 ) 2 + (l + B 2 ) -2] -2B In |b| -B(l + B 2 )T 

(i + b 2 ) 2L2|B| J ._ . 


where A and B remain as defined in equation (74). Thus, the integral 
which contains F in equation (70) is approximated as follows: 



Values of f x and g x are presented graphically in figures 1 and 2. 

(it is noted that the graphs in reference 21 that correspond to figure 1 
of this report are mislabeled.) 

With the simplifications introduced in this section, the function I 
of equation (56) is approximated by 


I(x) 




(79) 
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for Mach numbers less than those for which the shock wave is situated 
at the trailing edge of the airfoil, and 



for larger Mach numbers. 


Determination of 

The term up^ that appears in the integral equation for transonic 
flow represents the values for u given by linear theory for points on 
the airfoil surface. Its values can be obtained from the general two- 
dimensional solution for up given in equation (47) b y performing the 
indicated operations and setting z = 0. As noted in equation (5l), 
simplification occurs for the symmetrical nonlifting airfoils being con- 
sidered herein because 


Au = u u - Uj = 0 


thus 



81 P 3 8! p 3 d! 2 0 s 


(81) 
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o 


( 82 ) 



NACA TN 3096 


39 


The Cauchy principal value is understood in the last integral. The 
present calculations are for thin circular-arc airfoils for which the 
reduced ordinates of the upper surface are given by 



where t represents a reduced thickness ratio which is related to the 
actual thickness ratio as follows: 


r-i^i 

P 3 c 


m 


Performing the indicated operations gives the following expression 

v 



(85) 


We thus have, by substituting equation ( 85 ) into equation (59), 


L(x) - 2%, - 1 » f t 


1 + 



(86) 


It is seen that t plays the role of a similarity parameter. Thus, Ur 
for a family of symmetrical nonlifting airfoils having the same thick- 
ness distribution depends only on t and the position coordinate x/c« 
Inspection of the integral equations for transonic flow shows that their 
solutions for u^ also depend only on t and x/c. 

Many previous papers on transonic flow (e.g., refs. 12, 14, 16, 17, 
and 18) have used the symbol t Q to designate a Mach number thickness- 
ratio parameter different from the parameter f used herein. The defi- 
nition of ! q , suitably generalized to allow for various expressions 
for k, and the relation between £ Q and t are as follows: 


I 


o 


1-Mo 2 

[U Q k(t/c) ] 2/3 



(87) 
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The parameter 4 0 possesses the distinct practical advantage of 
approaching zero rather than infinity as the Mach number approaches 
unity. The corresponding parameter referring to the local conditions 
has commonly been designated by the symbol 4. 

4 = - ( 88 ) 

[Uok(t/c)] 2/3 

where M is the local Mach number. Since the quantity 1-M 2 is shown 
in equation (2lf) to be equivalent, in the present approximation, to 
1 -Mq -ku, the relation for I may be rewritten as 


1-Mq 2 - ku 

4 = 

[Uok(t/c)] 2/a 

which may be expressed in terms of u and f as follows: 


4 


u-1 


( 89 ) 


(90) 


Iteration Solution of Integral Equation 


Solutions of the simplified integral equations have been obtained 
for transonic flow about circular-arc airfoils by using a numerical and 
graphical process. Four slightly different techniques are used depending 
on the Mach number or, more precisely, on the value of t. One procedure 
is used for the subcritical Mach number range, a second, is used in an 
attempt to find supercritical shock-free solutions, a third to determine 
supercritical solutions in the range of f for which the shock position 
is forward of the trailing edge, and a fourth for still larger values 
of t. The latter range extends up to a free -stream Mach number of 
unity. 


Subcritical flows .- Solutions for small Mach numbers can be obtained 
by a direct iteration process starting with the linear-theory solution. 

In detail, the calculations proceed in the following manner. The values 
Of u L corresponding to a given t (specified thickness ratio and Mach 


number) are calculated from equation ( 85 ) . 


The 



curve is approxi- 


mated with a stepwise distribution of ten steps, and the values of the 
I curve are computed therefrom using equation (79). Knowing the values 
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of Uj^(x) and l(x), one obtains a first approximation to Uy(x) using 

equation (59)* The process is now repeated using the values for u^ 

to calculate a new I curve, from which a second approximation can be 

determined for uy, etc. A typical set of results illustrating the 

convergence of the process is shown in sketch (s). The process appears 

to converge rapidly, at least in the 

present calculations, whenever the 

Mach number (or more properly f) is 

sufficiently small that u does not 

exceed unity (or the local velocity 

does not become supersonic) at any 2 

point in the calculation. If t is u. 

increased to a value such that u 

becomes equal to one at any step in ' 

the iteration, however, the process 

terminates abruptly with the appear- 

ance of complex values of u. Since ? 

the results of successive iterations 
oscillate in this Mach number range 
and since the starting u^ distribu- 
tion provides maximum values that are Sketch (s) 

too small, any such termination occurs in the first iteration step. 

Such an abrupt termination of the calculations is in marked contrast to 
the well-known property of the classical iteration methods of producing 
second- and higher-order solutions indicating shock-free mixed flows. 

It will be shown in the following discussion that the difference in 
behavior is not a product of the additional assumptions and approxima- 
tions introduced herein, but stems directly from the recognition and 
retention of the quadratic nature of equation (6l). 

This point can be discussed in greater detail by confining atten- 
tion to only the first step of the iteration procedure in which the 
I curve is calculated using uy = u^. If the present iteration pro- 
cedure is used, if the u-^ curve is replaced with ten rectangular 

steps as shown in sketch (q), and if the I curve is calculated using 
equation (79 ) > then equation (66) yields the following result for Uy 
at the 50-percent-chord station of a circular-arc airfoil: 

*w(f) = l±/l - f T + 1 . 030 T 2 ( 91 ) 

For low Mach numbers, f is small, the discriminant is positive, and the 
first approximation for u w can be readily obtained using the minus 
sign. As t increases, the value of the discriminant decreases and 
becomes negative when t exceeds 0.490. This value of f, therefore, 
represents an upper limit for the Mach number at which the present 
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process can lead to a useful result if the starting Uy distribution 
is taken to be that given by linear theory. The value of u w (c/2) 
given by equation (91) is equal to ( Vf/rt) at small t , and attains a 
maximum value of unity when t = 0.490. 

The classical iteration methods of subsonic flow theory have not 
shown such an abrupt termination of solutions upon attainment of sonic 
velocity. The same situation occurs with the present equations if the 
iteration procedure is altered so as to conform with the classical 
methods. This change consists simply of substituting u w into not 

only the integral of equation (66) but also the term Uy 2 12 standing 
outside the integral, thereby converting the quadratic equation for uy 
into a linear equation. If this procedure is adopted together with the 
same value for the integral that was used to obtain equation (91), the 
result of the first iteration is 

* W (f) = n f +0 ' 295 f2 (92) 

Since a linear equation is solved at every step of the iteration process, 
the procedure never terminates. Values for uy. in the midchord region 
become larger with every iteration step, however, and it seems to be the 
prevailing belief that convergence is obtained only for Mach numbers less 
than the critical. 

It is of interest to compare the approximate results of the first 
iteration step reviewed above with the exact results for the same air- 
foil given by Hantzsche and Wendt in reference 32. If the latter results 
are made comparable to the present results by taking the limiting value 
corresponding to small-disturbance transonic flow theory, the exact 
result of the first iteration step is 

fi w (t)- ; ? + ( - j) f2 - ; f * 0-257 f 2 (93) 

The first term of both equations (92) and (93) is that given by linear 
theory and is the same in both calculations. The difference in the 
coefficient of the second term is the result of the. errors introduced 
in the approximate solution of equation (6l)(i.e., the velocity- 
distribution function, finite steps for u.y, etc.) and is some sort of 
a measure of the accuracy of the approximate calculations. Just as with 
equation (92), a value for uy can be calculated for all t, although 
the question remains of whether the result is a valid first step in a 
convergent process. 

If Hantzsche and Wendt had performed their iteration calculations 
in a manner comparable to that described herein so that the values of 
ul a re not introduced as an approximation for u in determining the 
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important influence of the points lying near P, they would have obtained 
the following relation (again to the approximation of transonic flow 
theory) as the result of the first iteration: 


5„(|) = 1 ± Jx - \ ?♦ (-! + 0f= . 1 ± fl - | f + 1 .: 


108 T 2 (94) 


This result is directly comparable with that' given in equation (91) 
except that the double integral of equation (6l) is evaluated exactly 
rather than approximately. Again, the accuracy of the approximations 
can be evaluated by comparing the two equations. More important, how- 
ever, is the fact that the exact results also terminate when t exceeds 
a certain value (O. 503 ) and that the corresponding value for u is 
unity. 


As noted above, the first step of 
the present iteration procedure start- 
ing with u L produces maximum values 
for u that are too large. Conse- 
quently, the calculations terminate at 
a value of f that is somewhat less 
than corresponds to the true critical 
Mach number. Solutions can be obtained 
for the remainder of the subcritical t 
Mach number range, however, if the 
starting uy is selected having the 
same general form but larger values 
than the corresponding u^ distribu- 
tion. A good starting distribution can 
be obtained by simply extrapolating the 
final results for smaller t. 

Calculations of the type just des- 
cribed have been carried through for 
several values of t less than the 
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. Sketch (t) 

value (tcj> = 0.598)> defined as the smallest value of f at 
which sonic velocity (uy = l) is attained in the final transonic solu- 
tion. The results of these calculations are presented in the form of 
ehordwise uy distributions for various t in sketch (t), and in the 


form of curves of uy versus t for various x/c in figure 3 / 


The 


-The results for f = f cr are also included on these graphs in order to 
complete the subcritical range, even though they are obtained using a 
^ a -^' eren ^ iteration procedure. Since the procedure is the same as that 
described in the next section, however, no further comment is necessary 
at this point. 
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corresponding values uj^ given by linear theory are also shown in 
figure 3. Comparison of the two sets of curves reveals that the values 
of u w coincide with those for u Lw at small t (low Mach number) but 
become larger than those for Ut at most stations on the airfoil as t 
approaches T cr . These trends, which are in accord with generally 
accepted experimental and theoretical results, will be discussed at 
greater length in a later section. 

Shock-free supercritical flows .- Although it is indicated in the 
Preceding discussion that the present iteration process does not lead 
to supercritical shock-free flow about airfoils when the starting point 
for the iteration calculation is u^* it is considered of importance 
and interest to ascertain the behavior of the present approximation when 
such a solution is ^deliberately sought by starting with sufficiently 
large values for u that I is greater than L. The interest in this 
matter is heightened by the fact that the less detailed method used by 
Oswatitsch leads to such results over a limited range of supercritical 
Mach numbers. (See figs. 4(c) and 4(d) of ref. 20 or figs. 7(c) and 7(d) 
of ref. 21.) Accordingly, iteration calculations have been performed 
starting with t greater than T cr and a symmetrical shock-free uy 
distribution containing a region of 'supersonic flow (u w > l) over the 
middle portion of the chord. A slight change in the iteration procedure 
is necessary, however, to eliminate the difficulty arising from the 
requirement described in an earlier section that I = L at the sonic 
points . To carry out the solution in the same manner as before would 
require that the uy distribution be found for which the corresponding 
I curve is tangent to a given L curve. Rather than attempting to 
find solutions by such indirect means, a method is used in which a new 
value of f is determined in such a manner that the tangency condition 
is fulfilled. In particular, the procedure is to calculate the I 
curve using equation (79) and the assumed values of t and uy. The next 
step consists of plotting the I curve and fitting an L curve computed 
from equation (86) for whatever value of t is necessary to fulfill the 
tangency condition as shown in sketch (j). A new set of values for u w 
can now be cal ciliated using equation (59). W 


uy = 1± a/i - L ( 59) 

where "the minus sign is used at stations upstream from the forward sonic 
point and downstream from the rear sonic point and the plus sign is used 
for the intermediate stations. 
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This process has been carried out for a number of assumed initial 
u w distributions. Sketch (u) shows a typical set of results obtained 



by starting with an initial value for t of O .65 and a u w distribu- . 
tion obtained by extrapolating the trends indicated by the solutions 
for subcritical flows. The assumed initial distribution is indicated 
by a dotted line in the upper portion of part (i) and the associated 
I and L curves are shown immediately below. The tangency requirement 
is fulfilled by taking f = 0.642 which compares well with the initial 
value of 0.65. The corresponding u w distribution calculated there- 
from is shown in the upper portion of part (i). Several points of 
interest are to be noted. The first, of course, is that the assumed 
and calculated distributions are distinctly dissimilar. The second 

is that the calculated distribution bears a marked resemblance to a 
result given by Oswatitsch (see fig. 4(d) of ref. 20 or fig. 7(d) of 
ref. 21) which may be described as a symmetrical u„ distribution con- 
taining an abrupt expansion on the front of the airfoil and a compres- 
sion shock symmetrically situated on the rear of the airfoil. 

Since the assumed and calculated u w distributions are too dis- 
similar to be regarded as solutions, the iteration process was continued, 
^si.ng the calculated results of the first step as the initial values for 
the second step. The results of this calculation are summarized in 
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part (ii). As indicated, t changes from 0.642 to 0.671 and the new uy 
distribution is again markedly dissimilar to the initial distribution. 

The calculations have been carried through several more steps and the 
results are shown in parts (iii) through (vi). It can be seen that the 
initial and calculated uy distributions do not agree, even qualita- 
tively, until several steps of the iteration process have been completed. 
In the meanwhile, the value of t has decreased to approximately that 
for the critical Mach number and the region of supersonic flow has effec- 
tively disappeared. The sixth and last step of the iteration calcula- 
tions presented in sketch (u) has produced a value of f = 0.600 and 
a u w distribution nearly identical with that shown in sketch (t) for 
the critical Mach number (? cr = 0.598). Although this sketch shows the 
results of only one series of calculations, similar results have been 
obtained starting with other assumed symmetrical, shock-free, super- 
c^*i^ical uy distributions. No case was found in which the calculated 
values repeated the assumed values until t had decreased to approxi- 
mately f cr and the supersonic zone had vanished. 

It should be noted before leaving this section that these results 
are not presented with the intention of proving or disproving anything 
about the more fundamental question of the existence or nonexistence of 
shock -free transonic flows. The purpose, rather, is to illustrate the 
behavior of the present approximation furnished by the simplified inte- 
gral equation and the iteration method of solution. 

Supercritical flows - shock wave forward of the trailing edge .- The 
preceding section has shown how the present method of calculating velocity 
distributions on thin airfoils fails to converge when t is greater than 
. ^cr and flow is assumed to be shock free. The identical iteration 
process will converge rapidly to a solution, however, if the initial u w 
distribution contains a discontinuity in accord with the shock relations. 
The method starts by selecting a value for f and assuming a reasonable 
distribution for Uy over the chord. The main point to observe in the 
selection of Uy is to include a shock wave (a discontinuity in uy) 
through which uy jumps from 1 + A immediately ahead of the shock to 
1 immediately behind the shock. As noted previously, such a jump in 
Uy is consistent with the assumption that the shock wave is a weak normal 
shock. The next step is to calculate the I curve using equation (79) 
and the assumed values of uy and f . The I curve is plotted and an L 
curve is computed using equation (86) and whatever value of t is nec- 
essary to fulfill the tangency condition, as illustrated in sketch (m) . 

A new set of values of u w can now be calculated using equation (59), 
taking proper care to change from the minus to the plus sign at the point 
of tangency and then return to the minus sign aft of the assumed shock 
position. In this way, a new approximation for the uy distribution 
corresponding to the new value for t is obtained, but the position of 
the shock wave is unaltered. If the new values for t and Uy are suffi- 
ciently close to the assumed values, it is presumed that an approximate 
solution has been found. In general, however, such a close correspondence 
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Sketch (v) 

is not obtained after the first step, 
and the entire calculation is repeated 
using the new values for f and uy 
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Sketch (w) 


in place of those as sumed initially. 

Throughout the process, the position 
of the shock wave is kept fixed, and 
the value of t is allowed to vary 
as necessary. In this way, the itera- 
tion process selects the solution for a given shock position rather 
than for a given Mach number. This process was carried out with the 
shock wave situated at several different stations on the airfoil chord 
and was found to converge rapidly even when the selected initial values 
for Uy and t differed considerably from their final values. A typical 
set of results (namely, that for the shock position fixed at 90-percent- 
chord station) is shown in sketch (v). In common with the rest of the 
calculations, the initial values for Uy and f were selected by extrap- 
olating the trends indicated by the calculations for more forward shock 
positions. 


Calculations of the type just described have been carried through 
with the position of the shock wave fixed at 60-, 70-, 80-, 90-, and 
100-percent chord. The results are presented as chordwise Uy dis- 
tributions for various t in sketch (w) and as the variation of uy 
with f at selected values of x/c in figure 4. 
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In order to test the method further , calculations were repeated in 
many of the cases using a different set of initial values for t and u w . 
The results of two such calculations for the shock position at the 
90 -percent-chord station are shown in sketch (x) . Although the results 




shown in part (i) were obtained using values for uy that were pur- 
posely selected to he too large and those of part (ii) to he too small, 
it can he seen that similar results are obtained after only a few steps 
of the iteration process. These same results are presented in a second 
manner in sketch (y) wherein the values of uy at several stations on 
the chord are plotted as a function of t. The final values to which 
the calculations converge (i.e., those given in fig. 4 ) are indicated by 
the solid line. The points connected by the dotted lines are the values 
obtained at each step of the calculations . It may be seen that the 
present procedure appears first to place the values of u w and t on the 
curve of correct solutions and, subsequently, to converge to the final 
result. 
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Supercritical flows - shock wave at the trailing edge . - The calcu- 
lations just described indicate that the shock wave moves rearward with 
increasing t and reaches the trail- 
ing edge when t = 1.34, correspond- 
ing to a Mach number of 0.92 for a 
4-percent-thick circular -arc section. 

The calculations for larger subsonic 
Mach numbers were performed under the 
assumption that a strong shock wave 
remains at the trailing edge. As in 
the preceding analysis, the shock 
wave is assumed to stand perpendic- 
ular to the x axis, and the flow is 
assumed to be normal to the shock wave. 

These latter conditions are not 
strictly correct at the base of the 
shock wave since the flow must turn 
through a finite angle. The error 
incurred, however, is small for thin 
airfoils . 



The method of calculation used 
for the higher Mach number range 
where the shock wave stands at the 
trailing edge is essentially the same 
as that used for mixed flows at lower 
Mach numbers. A slight modification 
enters in that the I curve is calcu- 
lated by means of equation (80) rather 
than (79) • This change is made because 

the contribution of the region imme- Sketch (y) 

diately behind the trailing edge becomes of increasing significance as 
the Mach number approaches unity. Since the contribution is small at 
the lowest Mach number for which the shock wave is situated at the 
trailing edge, no significant discontinuity is produced in the results 
by this change in procedure. A difficulty arises in the iteration 
process, however, because there is no longer any distinguishing feature 
to fix the value of f for which the solution is being sought. To 
review, for pure subsonic flows, f itself can be maintained at a fixed 
value from step to step in the iteration process. At Mach numbers some- 
what greater than the critical, where the shock wave stands on the air- 
foil, the tangency requirement makes it difficult to prescribe t directly 
but a satisfactory method is obtained by fixing the shock position and 
carrying out the iteration process until the associated value for t 
is found. At still larger Mach numbers, the first method is of no avail 
and the second method cannot produce a -unique result since the shock 
wave is considered to be fixed at the trailing edge over the rest of the 
range of subsonic free-stream Mach numbers. For lack of a better method. 
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the cal c\ilations were carried out for this range using a procedure that 
might he described as being partly iteration and partly trial and error. 
The calculations proceed as follows: A value of f is selected and an 

estimate made for the associated uy distribution. On the basis of 
these _values, calculations are performed resulting in a new set of values 
for t and uy. Similar calculations are repeated on a trial -and-error 
basis using the same Uy distribution but various values for t until 
the resulting t is equal to the assumed value. These values are then 
plotted on a graph of uy versus t. If these values were an actual 
solution of the integral equation, the resulting values of Uy would 
also be equal to the assumed values. In the present calculations, how- 
ever, the resulting values of u w are generally found to be somewhat 
smaller than the assumed values. These new values for u w , together 
with a smaller value for f, are next taken as the starting values for 
a second series of trial-and -error calculations . Again, the Uy distri- 
bution is held fixed as various values are tried for f. The calcula- 
tions are again terminated when the value of t is so selected that it 
repeats itself. The values of uy and t are then plotted on the graph. 
Unfortunately, it was not found possible in the present calculations to 
determine a set of values for uy and t that would repeat themselves 
precisely. In all cases investigated, the values of Uy and t were 
found to diminish somewhat in successive iteration steps. The source 
of this difficulty has not been ascertained at the present time. It 
could be due to one or more of the simplifying assumptions introduced 
in the approximate calculations or it could be due to the fact that an 
iteration process was not devised which would lead to a definite result, 
as at smaller f. 


There are a number of points , 
however, which tend to indicate that 
the successive values of uy and t 
obtained after the first few iteration 
steps may be regarded as near solu- 
tions. One of these concerns the fact 
that calculations made starting with 
different initial, values for uy and t 
converge to a common result after the 
first few iteration steps. To illus- 
trate, results of typical series of 
calculations starting with f of 
about 20 and three different uy dis- 
tributions are shown in sketch (z). 

It can be seen that the three sets of 
calculations all converge to determine 
a single line after the first few 
iterations. In order to promote 
insight into the significance of the 
indicated variation of t from step 
to step of the iteration process, the 
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corresponding Mach numbers for a 4-percent-thick section are also shown. 
Comparison on this basis shows that the variation of Mach number from 
step to step of the iteration process is very small. 

An additional factor lending credence to the applicability of the 
solutions for large f is provided by the phenomenon of the Mach number 
freeze wherein the local Mach number is invariant with changes in the 
free-stream Mach number when the latter is near unity or, more precisely. 



It is known from the papers of Vincenti and Wagoner (ref. 14), Liepmann 
and Bryson (refs. 28 and 29 ), and others that the corresponding approxi- 
mate relations yielded by the small -disturbance transonic theory is 



If the parameter k involved in the definitions of | is independent 
of Mo, as is the case when k is taken to be as given in equation (17), 
and if the local Mach number M is calculated using equation (24), the 
above two relations are completely equivalent. If k contains Mq, 
however, as it does in the preferred definition given in equation ( 7 ), 
dM/dMo only vanishes when M, as well as Mo, equals unity. If it is 
assumed that the freeze^ extends over a finite range of Mach numbers, the 
variation of u with t must obey the following relation for large f : 


i = 


u-l 

T 2 / 3 


= const. 


(97) 


Equation (97) has been used to compute a curve of u^ versus t under 
the assumption that the freeze extends at least to as small a value 
of t as 10, which corresponds to a Mach number of 0.978 for a 4-percent- 
thick section. These lines, which all have the characteristic form 


g W 1 

(SwW 1 



(98) 


are also shown in sketch (z). It can be seen that the lines so deter- 
mined are almost identical with those determined previously after the 
first few steps of the iteration process. This comparison shows that 
the solutions obtained by the iteration process possess the phenomenon 
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of the Mach number freeze. A further check was made by extrapolating 
the curves of u w versus t to t = 1+0 with the aid of equation ( 98 ) and 
using the values thus obtained as the initial values for an iteration 
calculation. It is found that these values nearly satisfy the integral 
equation although the values of u w and t diminish slightly from step 
to step during the iteration process, just as at the smaller values of f. 
The significance, in terms of Mach number, of these changes is even 
smaller than at lower Mach numbers. For example, the changes produced in 
the values of u and t in the first step of the iteration process corre- 
spond, for a 4 -percent -thick section, to an indiscernible change in the 
local Mach number and a change in the free-stream Mach number of less 
than 0.001. On the basis of these results, it appears proper to extrap- 
olate the curves of % versus t to values of f approaching infinity, 
corresponding to a Mach number of unity, by using equation ( 98 ). 



The results of the calculations 
for the range of f for which the shock 
wave is at the trailing edge are summa- 
rized in sketch (a') and figure 5 . 


RESULTS 

Results in Terms of Reduced Quantities 


The calculations described in the 
preceding sections have produced values 
for the velocity distributions at the 
surface of thin circular-arc airfoils 
in flows having free-stream Mach num- 
bers ranging from zero to unity. These 
results are presented in graphical form 
in figures 3 through 5 in terms of the 
reduced quantities u and t defined in 
equations ( 29 ) and (84) and repeated 
below 


“ **• _ 11 

U = U, T = 

1 ■ Mq 2 ; (1 - mq 2 ) 3/2 c 

where k represents the coefficient of 
the nonlinear term in the approximate 
differential equation for cp (eq. ( 6 )). 

The results of a number of pre- 
vious investigations of transonic -flow 
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theory have been presented in terms of the quantities ! 0 and i defined 
in equations ( 87 ) and (88) 

1 - Mo 2 1 - M 2 

= —> i = 

[U 0 k(t/c)] 2/3 [U 0 k(t/c)] 2 / 3 

The relations between the two sets of quantities are given by equa- 
tions ( 87 ) and ( 90 ) 


^o 



£ 


u— 1 
f 2/3 


Pressure distribution .- In many applications, the quantity of prime 
interest is the pressure distribution rather than the velocity or Mach 
number distribution. Because of the simple relation between the pres- 
sure coefficient and the perturbation velocity provided by equation (l6) , 
however, it is a simple matter to determine the pressure distribution 
once the velocity distribution is known; thus. 


2u = _ g(l-Mp 2 ) _ = _ 2(t/cf /3 / u \ 

U ° U o k U ( Uo k) 1/3 ^W 


( 99 ) 


The latter expression for Cp suggests the introduction of a reduced 
pressure coefficient Cp defined by 


C 


P 


«J 0 k > 1/3 

(t/c f' 3 


Cp - 



( 100 ) 


Equation (100) may be rewritten in terms of £ 0 and I as follows: 


C p = - 2(| - | 0 ) (101) 

The foregoing results have been used to calculate Cp for numerous 
stations on the airfoil chord for various (or f ) . The results so 

obtained are presented in figure 6 in the form of chordwise pressure 
distributions (C^ vs. x/c) for various £ 0 . This figure is presented in 
three parts: Part (a) contains the results for subcritical Mach numbers, 

part (b) for slightly supercritical Mach numbers for which the shock wave 
stands on the airfoil surface, and part (c) for still larger Mach numbers 
for which the shock wave is at the trailing edge. 
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The same results are presented in another manner in figure 7 in 
which the variation of Cp with £ 0 is plotted for various stations 
x/c along the airfoil chord. This form of presentation is the counter- 
part, in terms of reduced quantities, of the plots commonly found in 
many experimental investigations showing the variation of pressure coef- 
ficient Cp with free-stream Mach number Mq at prescribed points on 
the surface of wings and bodies. The corresponding curves given by 
linear theory are shown in figure 7 by dotted lines. These curves are 
computed using the equation 

obtained by direct substitution of equations (87) and (100) into equa- 
tion (85). It can be seen that the present results and those of linear 
theory are in good agreement for values of considerably less than 

zero. For £ Q near to zero, however, the present results display a 
behavior completely different from that indicated by linear theory. 

This is as it should be since it is well known that linear theory is 
totally inadequate for the analysis of steady two-dimensional flows 
about airfoils when the Mach number approaches unity, or, in the pres- 
ent notation, when £ Q approaches zero. 

In order to provide further ^information regarding the significance 
of the indicated variations of with £ Q , lines of constant £ have 
also been included in figure 7- The local velocities are subsonic if 
I is negative and are supersonic if £ is positive. The Mach number 
freeze is indicated in figure 7 by the curve representing the variation 
of Cp with at a given station x/c becoming parallel to a line of 
constant £. It can be seen from this figure that such behavior occurs 
at all stations on the airfoil for near zero. 

Pressure drag .- The foregoing paragraphs have been concerned .with 
the determination of pressure distributions on thin nonlifting circular- 
arc airfoil sections. Once these results are known, it is a simple 
matter to determine the section pressure drag coefficient c^ 


1 + 


{io2) 



(103) 


It is convenient at this point to introduce a reduced section drag coef- 
ficient eg defined in terms of reduced quantities, thus 



_ afzjf) 

Cp 


dx 


o 


die 


(104) 
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It is clear from the definitions of C p , t, z, and x that the relation 
between and cj is 


c d 


(u 0 k ) 1/3 2 r c „ , (a 0 k ) 1/3 


(105) 


It shoul d he remarked that some additional error is incurred in the 
present calculations of drag because the pressures are large and only 
poorly determined in the vicinity of the leading and trailing edges. 

The variation of c^ with ! 0 has been computed and the results 
are shown in sketch (b* ) . It can be seen that c^ is zero for i Q 
less than -1.4l8 (corresponding to the critical Mach number). The rapid 
rise of c^ as is increased beyond the critical is associated with 

the rearward movement of the shock and terminates abruptly when the shock 
reaches the trailing edge at ! Q = -0.825. The drag coefficient con- 
tinues to increase slowly with further increases of t Q although at a 



Sketch (b’) 

much reduced rate; and, finally, at a value of i Q somewhat less than 
zero, c£ becomes invariant with further changes in £ Q . This latter 
behavior is associated with the Mach number freeze, thus, 



d(Z u/ f ) 

dx = 

die 



(106) 
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Since the Mach number freeze corresponds, in the approximate theory to 
(d|/dl Q )t =0, equation ( 106 ) may be simplified to the folloving 

form: 0 


( dc &\ _ 4 P° d(Zti/f) 4 
\ai 0 4 o=0 ' C A dx 




= 4 


^h.E.” ^T.E. 


(107) 


where (Zn) T ™ and (Zvi) m „ refer to the ordinates of the upper sur- 

face at the leading edge and trailing edge, respectively. Since both 
of these quantities are zero for the circular-arc airfoils treated 
herein, it follows that 


£5|\ 

dl o4 0 =o 


= 0 


( 108 ) 


The corresponding results for symmetrical double-wedge profiles 
given in references 17, 13> and l4 by Trilling, Guderley and Yoshihara, 
and Vincenti and Wagoner, respectively, are also shown in sketch (b 1 ). 
Except for the value of c^ at = 0, the two sets of results bear 
only a qualitative resemblance. In some cases, the reason for the dif- 
ference is clear; in other cases, the reasons are more obscure. An 
example of the former concerns the result that the drag of the circular- 
arc section is zero for less than -1.408, whereas the drag of the 
double -wedge section remains finite for all £ 0 . The latter behavior 
results from the fact that the critical Mach number is zero for the 
double -wedge section. The reasons for the pronounced difference in 
shape between the two curves for drag at subsonic speeds are not so 
clear. The immediate explanation is that the shock wave moves rearward 
across the chord of the circular-arc airfoil at a more rapid rate than 
it does with the double -wedge airfoil. For the circular-arc airfoil, 
the shock wave has moved to the trailing edge when £ 0 equals -0.825* 
This condition marks the end of the rapid increase in drag. On the 
other hand, Trilling 1 s results for the double-wedge airfoil indicate 
that the shock wave does not reach the trailing edge until the free- 
stream Mach number is unity, or = 0. Whether this difference is an 
actual property of the solutions of transonic small -disturbance theory 
for these two profiles, or the result of simplifying assumptions intro- 
duced in either the present analysis or that of Trilling remains an 
unanswered question at the present time. 
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Results in Terras of Physical Quantities 


The preceding section has summarized the results of the present 
calculations of the pressure distribution and drag of thin circular-arc 
airfoils in flows having free -stream Mach numbers up to unity. These 
results are given in terms of the reduced parameters Cp, c^, and-! 0 > 
which possess the advantage of condensing the information for all thick- 
ness ratios onto a single curve, but the disadvantage of being somewhat 
complicated and unfamiliar. Consequently, it is the aim of this section 
to re-express these results in terms of the more conventional quanti- 
ties Cp, c^[, Mq, and t/c. 

Critical Mach number .- An airfoil property that is always of inter- 
est is the critical Mach number, M cr . The variation of the critical 
Mach number with thickness ratio can be readily determined from the 
result that £ cr = -1.408, thus, 


1 - M cr e 


=cr 


[U 0 k(t/c )] 2/3 


= -1.408 


(109) 


In this and in the remainder of the discussion, it is assumed whenever 
the results of numerical computations are presented that k is as 
defined in equation (7) and that y is /.OOt 
equal to 1.4. Sketch (c 1 ) shows a plot 
of the results of these calculations 
together with the corresponding results 
obtained using (a) linear theory and 
(b) thin-airfoil theory for incompres- 
sible flow together with the Karman- 
Tsien rule to account for the effect 
of compressibility. In both of the 
latter calculations, the variation 
of Cp with Mq was determined using 
equatiSn ( 23 ), as is customary in 
engineering practice. It can be seen 
that the present calculations indicate 
that the critical Mach number is some- 
what less than given by either linear 
theory or by the Karman-Tsien method. 

This is in accord with the results found Sketch (c 1 ) 

when the more exact theories are applied to thin airfoils. 



.06 .06 

t/c 


Pressure distribution .- The foregoing general results have been 
applied to determine the pressure distributions on a 4-percent-thick 
circular-arc profile. These results are presented in two alternative 
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forms in sketches (d' ) and (e 1 ). The first of these shows the chord- 
wise variation of Cp for various free-stream Mach numbers and illus- 
trates the development and rearward movement of the shock wave as the 



— Transonic theory 
- — Karman-Tsien 
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° M a 

Sketch (e T ) 


Sketch (d* ) 

free-stream Mach number increases 
beyond the critical. The second sketch 
shows the variation of Cp with Mg for 
selected points x/c on the airfoil 
— chord. Three sets of auxiliary lines 
\ x/c are also shown; they are (a) lines of 
\fO constant local Mach number computed 
using equations (24) and (25) and 
“ (b) lines showing the variation of Cp 
V \ B with Mq at selected points on the air- 
j foil chord computed using thin-airfoil 
q theory for incompressible flow together 
*\ ’ with the Karman-Tsien rule . At Mach 

numbers less than the critical, it can 
x 4 be seen that the pressure coefficient 
\ at a given point on the airfoil surface 
X 3 varies with Mach number in a manner 
v \ similar to that predicted by the Karman- 
\ Tsien rule. It is also apparent, as 
was pointed out earlier in connection 
r~ with equations (22) and (23) and 
n j sketch (b), that the value to be taken 
— -j for the critical Mach number differs 

slightly depending on whether the lines 
of constant local Mach number are com- 
puted by the present approximation 
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(eq. (24)) or by the exact isentropic relation (eq. (25)). The irregu- 
lar behavior of the pressures on the rear half of the airfoil at super- 
critical Mach numbers is associated with the passage of the shock wave. 
At Mach numbers near unity, a completely different behavior is evidenced 
in which the values of Cp change with Mq in such a manner that the 
local Mach number remains essentially constant. Direct calculation 
shows that the present theory implies the following approximate relation 
for dCp/dMQ at a free -stream Mach number of unity. 3 


( ^ 2 ) = Jt 2 

? +1 3 


(Cp) ' 
Mq=1 


k=H ° 2 (f) 


(no) 


This may be compared with the exact relation given by Vincenti and 
Wagoner in reference ik. 


(: 




dCp 


4 

7+1 


7+1 




(in) 


.024 


Pressure drag .- The variation of 
the section pressure drag coefficient 
c a with Mq has also been calculated 
for the percent -thick circular-arc 
airfoil and the results are shown in 1 
sketch (f f )« The general features of 
this curve are very similar to those 
discussed previously in connection with 
the corresponding curve of sketch (b ! ) 
for the reduced quantities. The major 
point of difference concerns the slope 
of the curve at a Mach number of unity i 
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Sketch (f) 


It might be noted here that numerous investigators have given the 
approximate relation 


\ dM Q y 


Mq=1 


7+1 ‘ 


k = 


7+1 


rather than that given in equation (llO). As indicated by the auxil- 
iary equations, the relations for dC p /dMo at Mq = 1 follow directly 
from the particular definition of k. It is interesting to observe 
that the expression for k given in equation (7) provides the more 
accurate relation for (dCp/dMo)^_ 1 , even though,- as noted in connec- 
tion with the discussion of equation (96), it provides the lesser 
accuracy in the determination of (dM/dM G ) M ^_ i . 
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Thus, the preceding discussion has disclosed that dc5/d§ 0 is zero at 
a Mach number of unity; whereas sketch (f 1 ) shows that dc^/dMo is 
negative at Mq = 1. Direct substitution gives the following value for 
this slope: 



( 112 ) 


Vincenti and Wagoner have shown in reference l4 that the exact relation 
for flow about a closed airfoil is 



(113) 


The negative value of the slope given in equation (112) arises from the 
fact that the quantity k which appears in all the reduced parameters 
(Cl, C£, i 0 , etc.) is a function of Mq. If, as in many other papers, k 
is equated to ( 7 + 1 ) /U 0 and is thus independent of Mq, the value of the 
slope is zero, thus 


■fa) = 0 , k = z±i 

' dM °^M 0 =l Uo 


(ilk) 


.Comparison With Experiment 


Inasmuch as the results described in the preceding sections were 
calculated after making numerous simplifications and approximations, 
not the least of which is the assumption of inviscid flow, it is desir- 
able to include some comparisons with experiment. This is particularly 
true for the present problem since it is well known that phenomena out- 
side the scope of potential theory, such as separation, boundary-layer 
shock -wave interaction, etc. are prominent features of transonic flow 
about airfoils. Since it is indicated in the preceding discussion that 
the present results are in general accord with the proven Prandtl-Glauert 
and Karman-Tsien results in the subcritical Mach number range, the 
following remarks will be confined to the supercritical range. 

There are at least three papers available which present results of 
detailed measurements of flow at high subsonic velocities about symmet- 
rical circular-arc airfoil sections, namely, reference 33 by Liepmann, 
reference 34 by Liepmann, Ashkenas, and Cole, and reference 35 by Wood 
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and Gooderum. The first two of these are concerned primarily with 
boundary-layer shock -wave interaction and contains statements casting 
doubt on the accuracy of the values indicated for the free -stream Mach 
number, a quantity of only secondary interest in their investigations. 
In general, it appears that values given for Mq in these two papers 
are somewhat too great. The more recent investigation of Wood and 
Gooderum appears to be better in this particular and possesses the 
advantage of being made with an interferometer so that knowledge is 
gained about the entire flow field. A disadvantage, however, is that 
the tests were made using a 12-percent-thick model, which severely 
strains the assumptions of transonic small-disturbance theory. Never- 
theless, the comparisons will be made with the data of reference 35* 

The studies of Liepmann (ref. 33 ) > Ackeret, Feldmann, and Rott 
(ref. 36 ), and others have shown that the boundary layer can have a 
profound influence on a transonic flow field. This immediately raises 
a question regarding the usefulness of a potential -flow theory, such as 
the present, which disregards the boundary layer completely. In order 
to illustrate better the nature of these effects, two interferograms 
of the flow about 12-percent-thick circular-arc airfoils are reproduced 
from reference 35 and shown in sketch (g 1 ) . The free-stream Mach 



laminar turbulent 

i n 

t/c = .12, M 0 - .88, R e s 600, OOO , ref. 35 

Sketch (g* ) 

number is 0.88, and the Reynolds number based on the chord is 600,000 
for both photographs. The conditions for the two flows differ in that 
the boundary layer is laminar in the flow pictured on the left and 
turbulent in that pictured on the right. The interferogram for the 
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laminar case shows that the shock waves are of the A type and that 
the flow separates near the midchord station. When the boundary layer 
is turbulent, however, it may he seen that the shock wave is of the 
simple single -wave type, and that the extent of the region of separated 
flow is greatly diminished. From this pair of photographs, it is appar- 
ent that two of the simplifying assumptions introduced in the course of 
the present analysis (i.e., (a) no flow separation occurs and (b) the 
shock wave is a single normal wave) are in better accord with the physi- 
cal phenomenon if the boundary layer is turbulent ahead of the shock 
wave than if it is laminar. Sketch (h 1 ) shows experimental pressure 




t/c=J2 , M 0 =.88 , R e = 600,000 , ref. 35. 

Sketch (h* ) 

distributions determined from the two interferograms of sketch (g*) 
together with the corresponding theoretical results. Because of diffi- 
culties in interpreting the interferograms, the experimental pressure 
distributions presented in reference 35 and reproduced in sketch (h l ) 
are terminated at the separation point. It can be seen that the theo- 
retical results are in substantial agreement with the experimental data 
available for the portion of the airfoil forward of the separation point. 
Although the experimental and theoretical pressure distributions matched 
equally well for all the test data of reference 35, the degree of corre- 
spondence must be attributed, in part, to a fortunate cancellation of 
errors since individual components involved in the perturbation analysis 
(e.g., the relation between w and dz/dx, Cp and the perturbation veloc- 
ity components, etc.) contain appreciable errors when applied to such a 
thick airfoil. Although the corresponding experimental values are not 
available for the pressures at stations aft of the separation point, it 
is presumed from other experimental data that the pressures reach their 
maximum negative value in the vicinity of the separation point and return 
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toward the free -stream value at stations farther aft on the airfoil, 
rather than varying in the manner indicated by the calculated results. 

Before leaving this topic, the following remarks should be made 
concerning experimental factors which may affect these conclusions. 

The first is that the authors of reference 35 do not consider the data 
for the turbulent boundary layer to be as reliable as that for the lami- 
nar boundary layer. The reason for this is that the method used to 
make the boundary layer turbulent (the upper half of the airfoil was 
mounted on a flat plate which extended one chord length forward of the 
leading edge of the airfoil) produced such a thick boundary layer that 
difficulties were encountered in correctly extrapolating the lines of 
the interferograms to the airfoil surface. 4 The differences in boundary 
layer thickness are clearly evident in the interferograms shown in 
sketch (g 1 )- The second stems from the fact . that the tests were con- 
ducted with a 12 -percent-thick airfoil. Comparisons of theory and 
experiment for such thick airfoils not only strain the small-disturbance 
assumptions of the theory but also emphasize unduly such features of 
the flow as the curvature of the shock wave and boundary-layer separa- 
tion which are disregarded completely in the theory. The third is con- 
cerned with the fact that results of a recent flight investigation by 
Harrin (ref. 37) have shown that, at Reynolds numbers of the order of 
20,000,000, there is very little difference in pressure distribution 
with laminar and turbulent boundary layers. In particular, it is found 
that the A type shock and large region of separated flow commonly 
associated with laminar boundary layers do not occur. If this result 
is substantiated by further experiments, it will be of particular sig- 
nificance in work such as the present where some assumption has to be 
made about the nature of the shock system, and the assumption that the 
shock wave is a single wave leads to the greatest simplification. 

The principal discrepancy between calculated pressure distributions 
and those measured with a turbulent boundary layer stems from the fact 
that the shock wave meets the airfoil surface at a station farther for- 
ward than calculated. The interferograms of sketch (g T ) show that a 
substantial portion of the forward shift is due to a pronounced curva- 
ture of the shock wave near the airfoil. This suggests a comparison of 
the calculated shock positions with not only the positions observed at 
the airfoil surface but also at some distance, say a half chord length, 
away from the airfoil. The results of such a comparison with the data 
for turbulent boundary layers given in reference 35 are shown in 

4 Results are also given in reference 35* although not for a Mach number 
of 0.88, in which the boundary layer is made turbulent by a wire 
stretched across the test section one chord length ahead of the lead- 
ing edge. When these results are plotted in the form shown in sketch 
(f 1 )* they fall about half way between the calculated results and 
those obtained with the half model mounted on a plate. 
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Theory Experiment sketch (i*). It can be seen that the 

tO\ r 1 7 — <3-, calculated shock positions vary with 

/ □ Mq in a manner which parallels the 

/ ’ { experimental results at z/c = 0.5 but 

/ /» , f not those at the airfoil surface. It 

9 ~f ” is almost certain that the details of 

/ this phenomenon are greatly influenced 

/ ^/’ry by separation and by interaction of the 

J~ , g - shock wave and boundary layer, and are, 

' hence, beyond the reach of separation - 
•§ / 11 Q free potential theory for transonic flow. 

^ / 0 s (7 <j> Since the pressure gradients and shock 

f strength increase with airfoil, thickness, 

| / Symbol z/c Source of turbulence it; mi 8 ht te presumed that these effects 

^ / 0 0 plate would be greater for thick airfoils 

C J □ ° 5 W 0 * fe than for thin airfoils. An indication 

/ & .5 wire of such a trend is furnished by the 

/ pressure-distribution data of Gothert 

/ (ref. 38) for NACA 0006, 0009, 0012, 

yg Q 2 Jo J Q 0015, and 0018 airfoils. Sketch (j') 

M 0 shows a plot of the variation with i Q 

. of the shock position at the surface of 

° e c ' each of these airfoils. If the shock 

positions were changing in accord with the similarity rules of transonic 
flow theory, these results would all fall on a single curve. It can be 
seen, however, that this is not the case and that the shock moves rear- 
ward across the chord at a slower rate for the thick airfoils than for 
the thin ones. It appears, therefore, that at best, the calculated shock 

positions will only agree with those 

1 1 1 1 — 7 1 1 1 1.0 found experimentally for very thin 

„/ airfoils. 


Sketch (i* ) 


Theory (circular- arc airfoil) 

-o— Experiment (NACA 0006 airfoil) 
- 0 - * * 0009 • 


Sketch (j*) 


_ g It is apparent that changes in the 

shock position will be accompanied, in 
many cases, by substantial changes in 
the pressure drag of the airfoil. In 
" 6 particular, the forward shift of the 
- shock noted for the thick airfoils will 
c diminish the region of high negative 
-4 pressures on the rear of the airfoil, 
thereby decreasing the drag. These 
effects will, of course, be zero until 
the critical Mach number is exceeded, 
will increase as the shock moves across 
the airfoil, and may diminish somewhat 
again as sonic free -stream velocity is 
jyO approached, since the actual, as well as 
calculated, shock position then approaches 
the trailing edge. These effects would 
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probably alter the curve for the pressure drag of a thick airfoil so 
that it would appear somewhat like , 

that illustrated in sketch (k 1 ). 


CONCLUDING REMARKS c d 


The foregoing results are 
encouraging in that they show that 
the introduction of a small number 
of rough, although reasonable. Sketch (k 1 ) 

assumptions leads to a relatively simple method for the calculation of 
pressure distributions on thin circular-arc airfoils at all Mach numbers 
up to unity. Perhaps the most important aspect of the present work is 
the recognition of the quadratic nature of the integral equation and the 
retention of this feature in the iteration solution. With the knowledge 
that acceptable results can be obtained without excessive effort, it 
appears worthwhile to re-examine the approximate solution of the equa- 
tions with an eye toward improvement, or elimination, of the simplifying 
assumptions. Probably one of the weakest elements of the present method 
is the velocity-distribution function introduced in equation ( 63 ). This 
particular function is used to determine the entire flow field but, 
actually, only insures that the velocity and velocity gradient have the 
correct value at the airfoil surface and that the velocity diminishes 
toward zero at infinity as l/z 2 . As can be seen by comparing the 
experimental and calculated results shown in sketch ( Z ’), this function 



Theory Experiment ref. 35 



Lines of constant C p , M 0 =.88, t/c=J2 

Sketch ( Z ' ) 



66 


NACA TN 3096 


succeeds reasonably well for circular-arc airfoils. On the other hand, 
such a simple function cannot be expected to give good results for all 
airfoils. For instance, the present velocity-distribution function 
cannot be expected to provide good results for airfoils having flat 
surfaces over a substantial part of the chord (e.g., wedge airfoils, 
etc.) since it indicates no attenuation with distance above an element 
of the airfoil where the radius of curvature of the surface in infinite. 
Gullstrand (ref. 22) has proposed a different velocity-distribution 
function which satisfies one more known requirement, but it has not been 
established as yet whether or not it is sufficiently general to cover 
all interesting cases. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Nov. 19 , 1953 
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APPENDIX 


PRINCIPAL SYMBOLS 


speed of sound 

speed of sound in the free stream 
critical speed of sound 

function defined in equations ( 63 ) and (64) 


pressure coefficient. 


(U 0 k) ] 


2 u ° 


(t/cf 


chord 


section pressure drag coefficient, 


(t/c)' 


M o p 

— U 0 2 c 


pressure drag 

function defined in equation ( 67 ) 
function defined in equation (7l) 
function defined in equations (72) and (73) 
function defined in equations ( 76 ) and (77) 
function defined in equations (57) and ( 58 ) 

coefficient of nonlinear term of differential equation for cp 
(See eqs . ( 7 ), (17), (l8), and (2l).) 


width of element used in approximating the chordwise velocity 
distribution 
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M 

M C r 

Mo 

n 

P 

P 

?o 

R 

Re 

r 2 

r 3 

s 

s 

t 

Uo 

u,v,w 

u 

V 

w 

IV (V ~ 

u,v,w 


local Mach number 

critical Mach number 

free -stream Mach number 

inward normal to surface enclosing R 

arbitrary point x,y,z 

static pressure 

free-stream static pressure 

region of integration 

Reynolds number 

J (x-i ) 2 + (z-n 2 
J (x - £) 2 + (y -fj) 2 + (z - l ) 2 

surface of shock wave 
wing semispan 

maximum thickness of profile 
free-stream velocity- 

perturbation velocity components parallel to x,y,z axes, 
respectively 



Cartesian velocity components on the two sides of a shock wave, 
with u being parallel to the flow direction ahead of the 
shock 
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X 


x,y,z 


S - x 
b 

Cartesian coordinates where x extends in the direction of the 
free -stream velocity 


x 

y 

z 

Z 

Z 


3 

7 

A 


^o 


I 

p c 

Z 

f 

<J> 

<P 


x 

Py 

Pz 

ordinates of wing profiles 

kU D 

Z 

P 3 

J 1 - M 0 2 

ratio of specific heats, for air 7 = 1. 14- 

difference between values of quantity on the upper and lower 
sides of the xy plane 

variables of integration corresponding to x,y,z 

_ 1-Mp 2 

[U 0 k(t/c )) 2/3 
1 - M 2 

[Uok(t/c )] 2 / 3 

free -stream density of air 
surface enclosing volume R 


velocity potential 
perturbation velocity potential 
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9 

♦ 


a 

b 

cr 

L 

L.E. 

I 

Mq = 

P 

S 

T.E. 

u 

W 


k_ 

P 2 


9 



Subscripts 

values ahead of shock wave 
values behind shock wave 

conditions associated with the critical Mach number 
values given by linear theory 
leading edge 

values on the lower surface of wing or wake 

values at Mq = 1 

values at arbitrary point P 

values along shock wave 

trailing edge 

values on the upper surface of wing or wake 
values at the wing surface 
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x/c 



T 

Figure 3.- Variation of u w and u Lw with t for various values of x/c. 

Suhcritical range. 



Figure 4.- Variation of u^ with t 
supercritical range, shock w; 


for various values of x/c 


ve forward of trailing edge 
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(a) t < 4.5 

Figure 5.- Variation of u w with f for various values of x/c. Medium 
and high supercritical range, shock wave at trailing edge. 
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(b) Supercritical - shock wave forward of trailing edge. 
Figure 6.- Continued. 
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Figure 6.- Concluded. 












